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CHAPTER  I 
INTRODUCTION 


1.0  OVERVIEW 

This  thesis  Is  based  upon  experimental  and  theoretical  work  In  the  area 
of  gear  mesh  failure  diagnostics.  More  specifically,  passive  diagnostic 
Instrumentation  was  installed  on  a  single  mesh  gear  fatigue  tester,  located  at 
NASA  Lewis  Research  Center,  to  periodically  record  the  gear  mesh  Induced  vibra¬ 
tion  v.om  Initiation  to  failure.  This  information  was  analyzed  using  several 
existing  gear  diagnostic  methodologies  to  determine  if  a  correlation  exists 
between  the  various  prediction  techniques  and  the  observed  modes  of  failure. 
This  thesis  presents  the  methods  used  and  results  obtained  when  applying  the 
diagnoctic  techniques  to  the  experimental  data: 

The  experimental  data  collected  during  this  thesis  consisted  of  discrete 
vibration  signatures  taken  from  the  eleven  gear  sets  that  were  run  to  failure. 
Two  accelerometers ,  both  with  a  frequency  range  of  0  to  10,000  Hertz,  and  an 
optical  sensor,  used  for  a  shaft  synchronous  signal,  were  Installed  on  an 
existing  spur  gear  fatigue  tester.  A  timer  was  Installed  so  that  the  vibration 
delta  could  be  periodically  recorded  on  an  FM  tape  recorder.  Of  the  eleven 
gear  sets  monitored,  five  failed  by  heavy  wear  and  scoring,  two  failed  by  sin¬ 
gle  tooth  pitting,  two  failed  by  tooth  breakage,  and  two  failed  by  distributed 
pitting. 

The  analytical  work  consisted  of  investigating  and  applying  several  gear 
mesh  failure  predication  techniques.  Two  of  the  methods  Investigated  were  the 
FM0  and  FM4  techniques  developed  by  Stewart  Cl).  FM0  is  a  general  method  used 
to  detect  a  variety  of  failures,  whereas  FM4  is  more  sensitive  to  a  single 
tooth  fault,  such  as  single  tooth  fatigue  cracks.  Also  used  was  the  technique 
utilizing  the  Hilbert  transform  to  demodulate  the  time  signal.  This  technique, 
proposed  by  McFadden  [23,  is  predominately  used  to  detect  fatigue  cracks  early 
In  their  development.  The  remaining  three  techniques  Investigated  were  the 
crest  factor,  sideband  level  ratio,  and  the  non-harmonic  to  harmonic  RMS  level 
energy  ratio.  [33 

This  thesis  is  divided  Into  five  chapters  to  help  separate  the  main  areas 
of  this  investigation.  The  remaining  two  sections  of  the  Introduction  present 
the  significance  of  condition  monitoring  In  rotorcraft,  and  a  brief  overview 
of  various  transmission  condition  monitoring  methods.  Chapter  II  covers  basic 
signal  processing  theory  and  techniques  used  In  this  analysis,  and  specific 
theory  on  each  of  the  diagnostic  methods  Investigated.  Chapter  III  addresses 
the  experimental  work  conducted  during  this  study  including  a  detailed  descrip¬ 
tion  of  the  various  gear  failures  found,  and  the  results  of  frequency  response 
measurements  performed  on  the  fatigue  test  rig.  Results  of  applying  the  vari¬ 
ous  gear  diagnostic  methods  to  the  data  collected  are  correlated  to  the  actual 
gear  failure  modes  In  Chapter  IV.  Finally,  Chapter  V  presents  the  conclusions 
of  this  research  and  recommendations  applicable  to  further  Investigations  Into 
gear  diagnostics. 

I. I  SIGNIFICANCE  OF  GEAR  CONDITION  MONITORING  IN  ROTORCRAFT 

In  aerospace  applications,  where  weight  and  size  are  premiums,  gear  sys¬ 
tems  are  usually  designed  close  to  their  projected  limits.  For  rotorcraft 
transmissions,  this  design  constraint  translates  into  frequent  transmission 
overhauls  and  frequent  transmission  related  accidents  resulting  In  death, 

Injury,  costly  damage,  and  fleet  grounding.  Current  on-board  condition  moni¬ 
toring  systems  not  only  provide  Insufficient  time  between  warning  and  failure, 


but  also  result  in  many  false  alarms  causing  unnecessary  and  costly  repairs. 
Thus  the  need  for  a  reliable  gear  train  condition  monitoring  system  is  para¬ 
mount  to  the  increased  safety  and  cost  efficient  operation  of  current  and 
future  rotorcraft. 

Some  specific  examples  can  be  given  illustrating  the  effects  of  transmis¬ 
sion  related  problems  on  fleet  readiness  and  safety.  The  entire  Marine  Corps 
fleet  of  93  CH-53E  Super  Stallion  Helicopters  were  grounded  In  February,  1987, 
due  to  a  defective  gear.  [4]  Also,  The  United  Kingdom  Department  of  Transport 
Air  Accidents  Investigation  Branch  < AAIB)  released  an  official  accident  report 
on  the  November  1986  accident  Involving  a  Boeing  234  Commercial  Chinook,  in  the 
North  Sea,  which  killed  45  people.  C5]  The  main  cause  of  the  accident  was 
determined  to  be  a  catastrophic  failure  of  a  modified  gear  In  the  forward 
transmission.  This  failure  led  to  the  desynchronization  of  the  twin  rotors 
such  that  the  forward  and  aft  rotors  clashed. 

The  early  detection  of  Incipient  failure  in  a  rotorcraft  transmission  is 
not  only  Important  for  safety,  but  also  overall  maintenance  Is  Improved  and 
becomes  more  cost  efficient.  Currently,  rotorcraft  transmissions  are  over¬ 
hauled  at  prescribed  intervals.  The  transmissions  are  overhauled  long  before 
the  calculated  life  of  any  of  the  critical  components.  Implementation  of  a 
reliable  on-board  condition  monitoring  system  can  Increase  the  time  between 
overhauls,  or  on-condltlon  maintenance,  without  sacrificing  safety  by  continu¬ 
ously  monitoring  the  state  and  wear  condition  of  critical  components.  The 
reduction  In  the  frequency  of  overhauls  not  only  decreases  maintenance  costs, 
but  also  decreases  the  probability  of  replacing  good  components  with  defective 
ones. 


1.2  OVERVIEW  OF  TRANSMISSION  CONDITION  MONITORING  METHODS 

There  are  currently  two  basic  methods  of  monitoring  the  condition  of  drive 
train  components.  The  first  uses  a  debris  monitoring  device  that  detects  the 
size  and  rate  of  wear  particles  in  the  transmission  lubricant  as  an  Indicator 
of  severe  wear  and  incipient  failure.  The  second  method  Is  based  upon  vibra¬ 
tion  data  obtained  using  one  or  more  sensors  mounted  on  the  transmission  case. 

Most  debris  monitoring  devices  use  magnetically  captured  debris  to  detect 
su  face  fatigue  failures  in  critical  gearbox  oil  wetted  components  such  as 
gears  and  bearings.  Some  of  the  debris  monitoring  devices  classify  the  cap¬ 
tured  particles  by  rough  sizes,  and  keep  account  of  the  rate  of  capture  and 
total  debris  count.  This  provides  an  Indication  as  to  the  damage  severity  of 
the  failing  component(s) .  One  problem  with  these  devices  Is  their  limited 
ability  to  detect  gear  failures,  since  in  most  cases  these-fai lures  do  not 
produce  large  amounts  of  metallic  debris  far  enough  In  advance  to  provide  suf¬ 
ficient  warning  time. 

The  use  of  vibration  analysis  methods  for  condition  monitoring  can  be  fur¬ 
ther  classified  Into  time  domain  and  frequency  domain  methods.  Time  domain 
methods  use  statistical  analysis  techniques  on  direct  or  filtered  time  signals 
to  detect  parametric  or  pattern  changes  as  transmission  components  wear. 
Statistical  methods  such  as  standard  deviation  and  kurtosis  are  used  to  qualify 
general  wear  from  tooth  specific  damage,  respectively.  Frequency  domain 
methods  use  the  Fast  Fourier  Transform  (FFT)  to  convert  the  time  signal  into 
its'  corresponding  frequency  components.  Vibration  energy  at  specific  frequen¬ 
cies  <1.e.  primary,  harmonics,  sidebands,  etc.)  can  be  used  to  monitor  gear- 
train  component  failures.  Both  methods  use  time  synchronous  averaging  to 
cancel  out  all  vibration  that  Is  non-periodic  with  the  shaft  frequency  being 
used  as  the  synchronous  signal. 
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CHAPTER  II 

THEORY  AND  TECHNIQUES 

An  understanding  of  the  theory  behind  the  diagnostic  techniques  Investi¬ 
gated  along  with  some  basic  digital  signal  processing  theory  was  needed  to  suc¬ 
cessfully  carry  out  this  analysis.  A  review  of  this  study  Is  presented  In 
this  chapter  to  give  the  reader  sufficient  background  Into  the  methods  used. 

2.0  DIGITAL  SIGNAL  PROCESSING  TECHNIQUES 

An  Important  part  of  the  experimental  portion  of  this  study  was  the  con¬ 
version  of  the  analog  signals  recorded  Into  digitized  data.  Digitizing  the 
data  allowed  the  techniques  being  Investigated  to  be  applied  to  the  data  effi¬ 
ciently  and  easily  on  a  standard  personal  computer.  Along  with  the  analog  to 
digital  conversion,  the  digital  Fourier  transform  was  used  to  convert  digitized 
time  data  to  frequency  data.  Because  of  the  significance  of  these  and  other 
digital  signal  processing  techniques  used  to  the  overall  analysis,  some  digital 
signal  analysis  theory  will  be  presented  below.  ! 

The  goal  of  the  analog  to  digital  conversion  process  is  to  obtain  this 
conversion  while  maintaining  sufficient  accuracy  In  frequency,  magnitude,  and 
phase  Information.  Two  basic  concepts  responsible  for  the  accuracy  of  the  con¬ 
version  are  sampling  and  quantization.  Sampling  Is  the  rate  of  time  signal 
digitization,  or  the  timing  between  the  digital  pieces  of  the  time  record. 
Quantization  Is  related  to  converting  the  analog  amplitude  value  to  a  digital 
value.  Of  the  two  concepts,  sampling  Is  usually  the  more  critical  because  It 
effects  frequency  as  well  as  amplitude  accuracies. 

All  analog  to  digital  (ADC)  conversion  devices  sample  at  constant  Incre¬ 
ments  of  time.  Two  theories,  Shannon's  Sampling  Theory  and  Rayleigh's  criteri¬ 
on,  relate  the  sampling  process  to  the  limits  of  the  Information  obtainable 
with  minimum  errors.  Shannon's  sampling  theory,  also  known  as  the  Nyqulst  cri¬ 
terion,  states  that  the  frequency  the  ADC  device  samples  must  be  greater  than 
two  times  the  maximum  frequency  of  Interest.  This  theory  outlines  the  highest 
frequency  detectable  In  a  signal,  with  minimum  errors,  given  a  specific  sam¬ 
pling  rate.  Raleigh's  criterion  is  related  to  the  ability  to  resolve  closely 
spaced  frequency  components.  For  a  given  time  record  of  length  T  (seconds), 
the  minimum  resolution  frequency,  or  the  lowest  frequency  component  measurable, 
Is  given  by  the  following  equation. 

f(mln)  .  |  (1) 

The  errors  associated  with  sampling  processes  can  be  grouped  into  two  cat¬ 
egories,  variance  and  bias.  Variance  Is  due  primarily  to  random  deviations  of 
the  sampled  signal  from  the  mean,  or  expected  value.  These  random  deviations 
are  often  caused  by  noise  In  the  signal  and  measurement  equipment.  The  easi¬ 
est  way  of  reducing  variance  errors  is  by  time  synchronous  averaging.  With 
this  method  the  signal  Is  averaged  In  the  time  domain  with  each  record  start¬ 
ing  at  the  same  point  In  a  cycle,  as  determined  by  a  synchronous  signal .  The 
portions  of  the  signal  not  coherent  with  the  synchronous  time  period,  such  as 
random  noise  and  signals  not  periodic  within  this  period,  are  averaged  out. 
Ideally  the  number  of  averages  Improves  the  signal  to  noise  ratio  by  the  fol¬ 
lowing  equation. 

S/N  (after  Nt  averages)  -  (fNt  >*<S/N)(or1glnal)  (2) 

where  :  S/N  •  signal  to  noise  ratio 
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Bias  errors  cannot  be  minimized  through  averaging  techniques.  Limitations 
of  the  equipment  or  techniques  used  are  the  prime  contributors  of  bias  errors. 
Aliasing,  a  common  bias  error,  results  from  the  limitation  of  using  finite 
Information  to  describe  an  infinite  process.  An  example  of  aliasing  can  be 
seen  In  In  Figure  2.1.  The  solid  line  plot  represents  the  actual  signal. 
Because  the  sampling  rate  was  conducted  at  a  frequency  lower  than  the  frequency 
of  Interest,  the  digital  representation  of  the  signal,  the  dashed  line  plot. 
Incorrectly  Implies  a  low  frequency  signal.  To  guard  against  aliasing  prob¬ 
lems,  Shannon's  sampling  theory  along  with  anti-aliasing  filters,  or  low  pass 
filters,  should  be  utilized  with  any  ADC  device.  The  anti-aliasing  filters 
serve  to  assure  that  any  frequencies  higher  than  the  frequency  of  Interest  are 
removed  from  the  signal  prior  to  the  ADC  process. 

One  of  the  most  powerful  tools  used  In  signal  analysis  is  the  Fourier 
transform.  The  basic  theory  behind  the  Fourier  transform  Is  that  any  continu¬ 
ous  time  signal  can  be  represented  by  a  series  of  sine  and  cosine  functions  of 
various  amplitudes,  frequencies,  and  phase  relationships.  Use  of  the  Fourier 
transform  produces  no  new  Information,  It  just  transforms  the  data  from  one 
Independent  variable  to  another.  In  most  Instances  time  domain  data1  Is  trans¬ 
formed  to  the  frequency  domain.  Determining  the  frequency  content  of  a  time 
based  signal  Is  especially  useful  when  analyzing  mechanical  systems  with 
rotating  components.  The  Fourier  transform  In  the  classical  case  Is  given  In 
the  following  equations. 


f<t)exp(-jwt)  dt  (forward  transform) 


F(w)exp(jwt)  dt  (Inverse  transform) 


where  :  f(t)«  function  in  the  time  domain 

F(w)«  function  In  the  frequency  domain 

Because  the  limits  of  Integration  of  the  classical  fora  of  the  Fourier 
transform  are  from  negative  to  positive  Infinity,  the  Disc  ste  Fourier  Trans- 
form  (DFT)  must  be  used  with  experimental  data.  Equations  fo.  tfco  DFT  are 
given  below. 


F(m).  1/N 


f(n)exp(-j2tron/N> 


(forward  transform) 


F(m)exp(j2»ran/N>  (Inverse  transform) 


The  OFT  Is  based  upon  a  set  of  assumptions  concerning  the  sequence  of  the 
discrete  data  values  collected.  These  assumptions  are:  1)  the  signal  must  be 
a  totally  observed  transient  within  the  time  period  of  observation,  or  2)  the 


signal  must  be  composed  only  of  harmonics  during  the  time  period  of  observa¬ 
tion.  Leakage,  a  bias  error,  Is  a  good  example  of  the  type  of  errors  resulting 
If  neither  of  the  assumptions  are  met  when  applying  the  DFT.  Figure  2.2  Illus¬ 
trates  the  resulting  transform,  labeled  "DISCRETE  FT",  as  compared  to  the  true 
transform,  labeled  "EXACT  FT".  The  time  signal  In  this  case  was  neither  tran¬ 
sient  nor  periodic  within  the  time  record,  thus  the  energy  In  essence  leaked 
out  into  other  frequency  lines  not  representative  of  the  actual  signal.  Win¬ 
dowing  Is  one  method  used  to  minimize  leakage  errors.  For  example,  the  Hanning 
window  multiplies  the  time  signal  by  a  function  that  gradually  reduces  to  zero 
at  the  end  of  the  time  period.  This  forces  the  signal  to  appear  to  be  tran¬ 
sient  within  the  time  period,  thus  satisfying  the  DFT  assumptions.  There  are 
functions  that  naturally  meet  the  DFT  assumptions,  these  are  called  self  win¬ 
dowing  functions.  These  functions  produce  no  leakage  errors  and  thus  require 
no  windowing. 

The  Fast  Fourier  Transform  (FFT)  was  developed  to  reduce  the  computational 
time  of  the  standard  DFT  process.  The  FFT  time  savings  results  from  partition¬ 
ing  the  DFT  matrix  Into  a  number  of  smaller  Individual  matrices  and  reordering 
the  data  for  more  efficient  computation.  This  results  In  the  computational 
time  required  for  the  FFT  to  be  a  function  of  Nlog2N  Instead  of  Ns  for  the 
DFT.  Along  with  the  time  saving  advantages  of  the  FFT  is  the  retention  of  the 
phase  Information,  which  permits  the  evaluation  of  multichannel  measurements 
such  as  transfer  functions,  and  coherence.  [6] 

The  frequency  response  function  Is  another  Important  digital  signal  pro¬ 
cessing  method  which  Is  used  to  describe  the  dynamic  properties  of  a  physical 
system,  The  classical  definition  of  the  frequency  response  function  is  the 
Fourier  transform  of  the  unit  Impulse  response  function,  as  shown  in  the  fol¬ 
lowing  equation.  [7] 


H(f>- 


h<t>expH2*ft>  dt 


(5) 


where  :  H(f>«  Frequency  response  function 

h<t>»  Unit  Impulse  response  function 


The  frequency  response  function  gives  a  direct  relationship  of  the  output  to 
the  input  of  a  system  In  the  frequency  domain.  This  relationship  Is  especially 
useful  when  trying  to  determine  how  a  signal  Is  altered  as  It  passes  through  a 
system. 

There  are  several  ways  of  estimating  the  frequency  response  function  using 
measured  quantities.  The  first  method  Is  simply  computed  by  dividing  the 
Fourier  transform  of  the  output  by  the  Fourier  transform  of  the  input,  as  given 
in  Equation  6. 


Hyx(f>»  Y(f)/X(f)  (6) 

where  :  Hyx  estimates  H 

Y(f )  and  X(f >  are  the  Fourier  transforms 
of  the  output  and  Input,  respectively. 
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To  reduce  the  variance  error  In  the  estimate  of  the  frequency  response  a  number 
of  samples  must  be  Involved,  as  seen  In  Equation  7. 


Y<f ), 


Hyx(f) 


<7> 


X(f>, 


The  above  estimates  are  reasonable,  however,  to  obtain  the  optimum  esti¬ 
mate  of  the  frequency  response  function  In  the  presence  of  noisy  signals,  the 
least  squares  technique  Is  often  used.  This  technique  produces  an  unbiased 
estimate  of  the  frequency  response  function  even  In  the  presence  of  Input 
noise,  assuming  that  the  Input  signal  and  the  Input  noise  are  not  correlated. 
[7]  The  frequency  response  estimate  In  this  case  Is  given  by  Equation  8. 


Hyx(f)  .  Gyx(f)  /  Gxx(f) 


(8) 


where:  Gyx(f) 


NS 


1.1 


Y(f)|  X*(f> 


Gxx(f)  - 


NS 


1-1 


XCf),  X*<f) 


1 


In  this  quatlon  Gyx  represents  the  cross  power  spectrum,  and  Gxx  represents 
the  spectrum  of  the  Input,  both  of  which  can  be  calculated  even  In 

the  presence  of  noise.  The  phase  Information  Is  preserved  In  the  term  Gyx. 

The  cross  and  auto  power  spectrums  are  also  used  to  calculate  another  Important 
signal  processing  parameter,  the  coherence  function. 

The  coherence  function  provides  a  means  of  determining  the  causality  rela¬ 
tionship  between  the  input  and  output  of  a  system.  The  coherence  function  Is 
given  In  Equation  9  below. 


where: 


v2(f> 


Gxy<f> 

GxyMf) 

Gxx(f) 

Gyy(f) 


Gxv(f)  Gxv*(f) 
Gxxiff  >  Gyy<f> 1 


(9) 


The  cross  spectrum 

The  complex  conjugate  of  the  cross  spectrum 
The  auto  power  spectrum  of  the  Input 
The  auto  power  spectrum  of  the  output 


The  coherence  function  is  a  measure  of  the  amount  of  the  output  signal  that  1$ 
directly  related  to  the  Input  signal  at  any  specific  frequency.  If,  for 
example,  at  a  certain  frequency  the  coherence  function  has  a  value  of  1.0,  then 
the  output  Is  a  direct  result  of  the  Input  signal  at  that  frequency.  If  the 
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coherence  function  has  a  value  less  than  unity  but  greater  than  zero,  then  one 
or  more  of  the  following  conditions  may  exist,  as  outlined  by  Bendat  [7]: 

1)  Extraneous  noise  Is  present  In  the  measurements. 

2)  Resolution  bias  errors  are  present  In  the  spectral  estimates. 

3)  The  system  relating  y(t)  to  x(t)  Is  not  linear. 

4)  The  output  y(t>  Is  due  to  other  Inputs  besides  x(t). 

The  coherence  function  Is  normally  used  In  conjunction  with  frequency  response 
measurements.  Here  the  coherence  function  provides  a  measure  of  the  quality  of 
the  measurement  by  giving  the  relationship  of  the  output  signal  to  the  Input 
Impulse.  Another  way  of  representing  the  coherence  function  Information  Is  by 
the  signal  to  noise  ratio.  The  signal  to  noise  ratio  Is  easily  found  using  the 
coherence  function  as  shown  In  Equation  10. 

2 

S/N  (f)  .  Y  - -  (10) 

1  -  Y*(f> 


2.1  THEORY  OF  GEAR  FAILURE  PREDICTION  METHODS 


Several  gear  failure  prediction  methods  were  Investigated  and  applied  to 
the  experimental  data.  The  basic  theory  behind  each  method  Is  given  In  the 
following  sections. 


FMO  METHOD 

The  FMO  parameter  Is  a  time  domain  discriminant  proposed  by  Stewart  Cl 3 
that  provides  a  simple  method  to  detect  major  changes  In  the  meshing  pattern. 
FMO  is  defined  as  the  ratio  of  the  peak-to-peak  level  of  the  signal  average  to 
the  sum  of  the  rms  levels  of  the  meshing  frequency  and  Its  harmonics,  as  given 
In  Equation  1 1 . 


FMO  - 


<U> 


where:  PP  -  peak-to-peak  level  of  signal  average 
A  -  amplitude  at  mesh  frequency  (1-1)  and 
harmonics  <1>1) 


FMO  Is  formulated  to  be  a  robust  Indicator  of  major  faults  In  a  gear 
mesh.  Major  tooth  faults,  such  as  breakage,  usually  results  In  an  Increase  In 
peak-to-peak  level  and  no  appreciable  change  In  meshing  frequency  levels, 
which  causes  FMO  to  increase.  For  heavy  uniform  gear  wear,  the  peak  to  peak 
level  remains  somewhat  constant  and  the  meshing  frequency  levels  tend  to 
decrease,  causing  the  FMO  parameter  to  Increase.  In  heavy  wear  situations, 
the  gear  tooth  profile  and  surface  quality  degrades,  causing  the  meshing  energy 
to  be  redistributed  from  the  meshing  frequencies  to  the  modulating  sidebands. 
From  the  physical  viewpoint,  two  surfaces  In  sliding  contact  produce  more  ran¬ 
dom  fluctuations  In  the  signal  as  the  amount  and  magnitude  of  their  surface 
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variations  Increase.  This  explains  the  redistribution  of  energy  In  the  fre¬ 
quency  domain  with  the  Increasing  wear.  For  minor  tooth  damage,  neither  peak- 
to-pe?.k  levels  or  meshing  frequencies  levels  change  an  appreciable  amount, 
thus  FMO  response  In  this  case  Is  limited.  E8] 

FM4  METHOD 

FM4  was  developed  to  detect  changes  In  the  vibration  pattern  resulting 
from  damage  to  a  single  tooth.  FM4  analysis,  proposed  by  Stewart  [1],  filters 
out  the  regular  meshing  components  from  the  signal  average  and  performs  two 
statistical  operations,  standard  deviation  and  kurtosls,  on  the  difference  sig¬ 
nal.  Equation  12  and  Figures  2.3  through  2.5  Illustrate  the  formulation  of  the 
difference  signal . 

D(t)  -  A(t>  -  R(t)  (12) 

where:  A(t)  »  original  signal  average 

R(t>  ■  regular  meshing  components  of  signal  average 
D(t)  •  difference  signal 

Figure  2.3  represents  a  plot  of  an  actual  signal  A(t)  after  time  synchro¬ 
nous  averaging.  Figure  2.4  represents  a  plot  of  the  regular  meshing  components 
R(t)  of  that  signal.  R(t)  Is  found  by  taking  the  FFT  of  the  original  signal, 
extracting  the  regular  components  and  taking  the  Inverse  FFT  of  these  compo¬ 
nents.  The  regular  components  consist  of  the  shaft  frequency  and  Its  harmonics, 
the  primary  meshing  frequency  and  Its  harmonics  along  with  their  first  order 
sidebands.  The  first  order  sidebands  are  generally  due  to  runout  of  the  gear 
because  of  machining  or  assembly  Inaccuracies  C93,  and  thus  are  considered 
regular  components.  The  difference  signal  Is  then  found  by  subtracting  the 
regular  components  from  the  original  signal.  Figure  2.5  shows  a  plot  of  the 
resulting  difference  signal.  The  FM4  analysis  method  Is  comprised  of  the 
standard  deviation  and  kurtosls  of  the  difference  signal  along  with  a  squared 
representation  of  the  difference  signal.  The  square  of  the  difference  signal 
represented  In  Figure  2.5  Is  shown  In  figure  2.6.  As  seen  In  this  figure  the 
square  of  the  difference  signal  magnifies  any  abnormalities  present  In  the  dif¬ 
ference  signal. 

The  FM4  analysis  method  uses  standard  deviation  and  kurtosls  to  extract 
Information  from  the  resulting  difference  signal .  The  standard  deviation  of 
the  difference  signal  Indicates  the  amount  of  energy  In  the  non-meshing  compo¬ 
nents,  where  the  kurtosls  Indicates  the  presence  of  peaks  in  the  difference 
signal.  The  theory  behind  FM4  Is  that  for  a  gear  In  good  condition  the  dif¬ 
ference  signal  would  be  primarily  noise  with  a  Gaussian  amplitude  dlstr'butlon 
[83.  The  standard  deviation  should  be  relatively  constant,  and  a  normalized 
kurtosls  value  of  three  would  be  reflected  [33.  When  a  tooth  develops  a  major 
defect,  such  as  root  cracks  or  severe  spalling,  a  peak  or  series  of  peaks 
appear  in  the  difference  signal.  The  kurtosls  value  would  thus  increase  to 
reflect  these  peaks.  The  standard  deviation  will  increase  only  when  the 
peak(s)  become  severe  enough  to  bring  up  the  rras  level  of  the  entire  difference 
signal.  It  should  be  noted  that  the  standard  deviation  will  also  Increase  with 
Increases  In  uniform  gear  wear.  This  Is  to  be  expected,  since  the  standard 
deviation  of  the  difference  signal  is  an  indicator  of  the  energy  level  of  the 
non-meshing  components,  which  Increases  as  the  gear  wears.  Thus  both  the 
standard  deviation  and  the  normalized  kurtosls  should  be  used  In  conjunction 
with  each  other  to  detect  the  onset  and  severity  of  single  tooth  defects. 


Standard  deviation  Is  defined  as  the  second  statistical  moment  of  an 
array  of  values  about  the  mean  of  those  values.  It  Is  an  Indicator  of  the  var 
lance  of  the  values  In  the  array.  The  standard  deviation  of  the  difference 
signal  can  be  found  using  the  following  equation. 


RMSDS 


1/2 


(13) 


where  :  RMSDS  -  Standard  Deviation 

d  -  Mean  value  of  the  signal 


Kurtosls  Is  defined  as  the  fourth  statistical  moment  of  an  array  of  values 
about  the  mean  of  those  values.  It  Is  an  Indicator  of  the  existence  of  major 
peaks  In  the  array.  The  digital  form  of  the  kurtosls  equation  Is  given  In 
Equation  14. 


K  rn  1/N 


(14) 


where:  K«  Kurtosls 

d  -  Mean  value  of  signal 

This  absolute  kurtosls  value  will  Increase  proportionally  with  general  increases 
in  the  standard  deviation.  To  keep  the  kurtosls  parameter  sensitive  to  single 
tooth  damage  only,  the  normalized  kurtosls  equation  Is  given  below. 


US) 


where:  NK  ■  Normalized  kurtosls 
d  *  Mean  value  of  signal 


As  seen  In  Equation  15.  the  normalized  kurtosls  Is  simply  the  absolute 
kurtosls  divided  by  the  fourth  power  of  the  standard  deviation,  or  the  square 
of  the  variance.  This  allows  the  normalized  kurtosls  to  be  sensitive  only  to 
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peaks  In  the  difference  signal,  regardless  of  changes  In  the  standard 
deviation . 

Because  the  normalized  kurtosls  values  are  non-dimensional*  signals  with 
different  magnitudes  but  similar  shapes  will  have  similar  values.  A  square 
wave  is  found  to  have  a  normalized  kurtosls  value  of  1.0,  for  a  sine  wave  the 
value  is  1.5,  and  for  a  signal  of  essentially  noise  with  a  Gaussian  amplitude 
distribution  the  normalized  kurtosls  value  Is  found  to  be  3.0  [3].  Thus,  a 
normalized  kurtosis  value  greater  than  3.0  Is  Indicative  of  a  peak  or  series  of 
peaks  existing  In  the  signal.  One  pitfall  with  the  normalized  kurtosls  parame¬ 
ter  Is  Its  drastic  decrease  in  peak  sensitivity  as  the  number  of  peaks  of  simi¬ 
lar  magnitudes  increase  beyond  two.  Thus  for  failures  Involving  two  or  more 
teeth,  the  normalized  kurtosls  value  may  not  Increase  far  beyond  3.0,  and  the 
failure  must  be  detected  with  the  standard  deviation  level. 

HILBERT  TRANSFORM  METHOD 

A  technique  was  developed  by  McFadden  [2]  to  detect  local  gear  defects, 
such  as  fatigue  cracks,  using  the  Hilbert  transform.  The  basic  theory  behind 
this  technique  Is  that  the  sidebands  around  the  dominant  meshing  frequency 
modulate  the  meshing  frequency  to  produce  the  time  average  signal.  Using  the 
Hilbert  transform,  the  signal  can  be  demodulated,  resulting  In  the  correspond¬ 
ing  amplitude  and  phase  modulation  functions.  The  phase  modulation  function 
is  especially  sensitive  to  fatigue  cracks  by  Indicating  a  phase  lag  at  the 
point  the  cracked  tooth  goes  into  mesh  [23. 

The  Hilbert  transform  is  primarily  used  to  transform  a  real  time  signal 
into  a  complex  time  signal  with  real  and  Imaginary  parts.  The  real  part  of 
the  complex  time  signal  Is  the  actual  time  signal,  and  the  Imaginary  part  Is 
the  Hilbert  transform  of  the  actual  time  signal.  This  complex  time  signal  Is 
referred  to  as  the  analytic  signal,  and  Is  given  in  the  following  equation. 


AN(t>  -  A(t>  ♦  tHCA(t) 1  (16) 

where:  AN(t)  -  analytic  signal 

A(t)  *  original  signal 

H ( A < t > 3  *  Hilbert  transform  of  origintl  signal 


The  Hilbert  transform  of  a  real  valued  time  signal,  A(t),  is  defined  as  the 
convolution  of  the  time  signal  with  l/*t,  as  shown  In  Equation  17. 


H[A(i>] 


A<T>nm-t))  dx 


(17) 


An  easier  way  of  computing  the  Hilbert  transform  is  to  utilize  the  convolution 
theorem,  i.e.  the  convolution  of  two  signals  in  the  time  domain  Is  equivalent 
to  the  inverse  Fourier  transform  of  the  product  of  the  two  signals  in  the  fre¬ 
quency  domain.  Thus,  to  determine  the  Hilbert  transform  of  a  real-valued  time 
signal  the  signal  must  be  transformed  to  the  frequency  domain,  phase  shifted 
by  -90  degrees,  and  transformed  back  to  the  time  domain.  The  determination  of 
the  Hilbert  transform  using  this  method  Is  illustrated  in  Equation  18. 
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(18) 


H[A(t)3  =  F-Uc-1  sgnf)  A(f)] 

where:  A(f)  =  Fourier  transform  of  original  signal 
sgnf  =  1  for  f>0  ,  -1  for  f<0 
F*k]  »  Inverse  Fourier  transform 


Once  the  analytic  signal  Is  found,  the  amplitude  and  phase  modulation 
functions  can  be  determined-  The  amplitude  modulation  function,  or  envelope 
of  the  signal,  Is  the  magnitude  of  the  analytic  signal,  and  is  given  In 
Equation  19. 


|AN(t) |  *  [(A(t)>2  +  (H[A(t>3>231/2  (19) 


The  phase  modulation  function,  or  Instantaneous  phase  variation.  Is  found  by 
using  Equation  20  below. 

<j>(t)  »  tan-1CHCA(t)]/A(t)]  -2irf0t  (20) 

where:  f0  *  carrier  frequency 


The  phase  modulation  function,  $(t),  represents  the  Instantaneous  phase  angle 
variation  with  respect  to  the  nominal  gear  speed  [10].  The  second  term  In 
Equation  20  represents  a  ramp  function  with  a  frequency  equal  to  the  carrier 
frequency,  f0,  that  Is  being  modulated.  This  term  Is  required  to  separate 
the  Instantaneous  phase  angle  variations  from  the  constant  carrier  frequency 
phase  function. 

Before  creating  the  analytic  signal,  the  original  signal  must  be  bandpass 
filtered  about  a  dominant  meshing  frequency.  This  dominant  frequency  Is  either 
-the  primary  mesh  frequency  or  one  of  Its  harmonics,  whichever  appears  to  give 
the  most  robust  group  of  sidebands.  The  width  of  the  bandpass  filter  depends 
on  the  location  of  the  meshing  frequency  to  other  meshing  frequency  harmonics. 
McFadden  [23  suggests  using  a  bandwidth  giving  the  maximum  amount  of  sidebands, 
even  If  they  Interfere  with  the  sidebands  from  other  harmonics.  The  reasoning 
here  is  to  Include  as  much  of  the  primary  modulating  sidebands  as  possible, 
assuming  that  the  Interference  from  the  other  sidebands  Is  negligible.  The 
dominant  frequency,  or  carrier  frequency,  Is  then  removed  and  the  resulting 
bandpassed  sidebands  are  used  In  constructing  the  analytic  signal  and  modula¬ 
tion  functions. 

As  stated  earlier,  the  phase  modulation  function  appears  to  be  the  most 
sensitive  of  the  two  modulating  functions  to  gear  fatigue  cracks.  One  particu¬ 
lar  feature  of  the  phase  modulation  function  Is  Its  ability  to  distinguish 
between  an  actual  fatigue  crack  and  a  particle  on  the  tooth.  For  a  fatigue 
crack,  the  phase  modulation  function  would  experience  a  phase  lag  at  the  point 
the  cracked  tooth  went  Into  mesh.  For  a  particle  on  the  tooth,  the  phase  modu¬ 
lation  function  would  experience  a  phase  lead  when  that  tooth  went  Into  mesh 
[23. 
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CREST  FACTOR 

The  crest  factor,  CF,  Is  a  simple  measure  of  detecting  changes  to  the  sig¬ 
nal  pattern  due  to  Impulsive  vibration  sources,  such  as  tooth  breakage  C33 . 

The  crest  factor  Is  easily  calculated  by  dividing  the  peak,  level  of  the  signal 
average  to  the  standard  deviation  of  the  signal  average,  as  given  In 
Equation  21 . 


r,  PL  (21) 

w  *  RMS 

where:  PL  «  peak  level  of  signal 

RMS  »  standard  deviation  of  signal  average 


Peaks  In  the  signal  will  result  In  an  Increase  In  CF.  . 


SIDEBAND  LEVEL  FACTOR 


The  sideband  level  factor,  SLF,  Is  a  course  indicator  of  single  tooth 
damage  or  gear  shaft  damage.  [81  To  calculate  SLF  the  first  order  sideband 
levels  about  the  primary  meshing  frequency  are  divided  by. the  standard  devia¬ 
tion  of  the  signal  average,  as  seen  in  Equation  22. 


SLF 


where:  FOSL  -  first  order  sidebands  level 
RMs  ■  standard  deviation  of  signal 


(22) 


A  bent  or  damaged  shaft  will  result  in  an  eccentric  meshing  pattcfn,  directly  \ 
reflecting  an  Increase  In  the  first  Order  sidebands  level >  thus  increasing  the 
SLF  value.  It  Is  unclear  to  the  author  how  this  factor  is  sensitive  to  single 
tooth  damage. 


ENERGY  RATIO 

The  energy  ratio  <£R)  is  formulated  to  be  a  robust  Indicator  of  heavy  uni¬ 
form  wear.  This  factor  divides; the  standard  deviation  of  the  difference  signal 
by  tne  standard  deviation  of  the  signal  composed  of  the  regular-  frequency  com¬ 
ponents  only  C31.  The  difference  signal  and  the  regular  components  signal  are 
the  same  as  those  defined  previously  la  the  FM4  analysts,  section.  Equation  23 
illustrates  the  ER  factor. 


ro  RH5PS 

tK  “  RMSSC 


(23) 


where;  RM|DS  »  standard  deviation  of  difference  signal  > 
RMSRC  «  standard  deviation  of  regular  components 
:  portion  of  the  signal 
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For  heavy  wear,  the  mesh  energy  redistributes  from  the  regular  meshing  fre¬ 
quency  components  to  the  modulating  sidebands.  This  would  result  In  an 
Increase  In  the  rms  level  of  the  difference  signal,  the  numerator,  and  a 
decrease  In  the  rms  level  of  the  regular  components  signal,  the  denominator. 
The  net  result  would  be  an  Increase  In  the  ER  value  In  a  more  robust  fashion 
than  just  the  Increase  In  the  difference  signal  rms  level,  as  calculated  In 
the  FM4  analysis. 


Amplitude 
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Figure  2.5  :  Plot  of  Difference  Signal,  D(t),  for  One  Shaft 
Revolution. 
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Figure  2.6  :  Plot  of  Squared  Difference  Signal,  [D(t)j^, 
for  One  Shaft  Revolution 
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chapter  III 

experimental’  procedures 


The  experimental  portion  of  this  study  was  conducted  on  a  gear  fatigue 
testing  apparatus  at  NASA  Lewis  Research  Center.  Some  background  Information 
on  the  fatigue  tester  will  be  presented  In  this  chapter  along  with  some  details 
on  the  type  of  gears  tested  and  their  specific  modes  of  failure.  Also  dis¬ 
cussed  In  this  chapter  Is  the  test  Instrumentation  used,  basic  test  set-up,  and 
some  frequency  response  measurements  performed  on  the  fatigue  tester. 

3.0  SPUR  GEAR  FATIGUE  APPARATUS 

The  spur  gear  fatigue  test  apparatus  at  NASA  Lewis  Research  Center  first 
became  operational  In  1972.  It  has  supplied  the  means  to  obtain  crucial  data 
on  the  effects  of  gear  materials,  gear  surface  treatments,  and  lubrication 
types  and  methods  on  the  fatigue  strength  of  aircraft  quality  gears.  A  cut 
away  view  and  schematic  diagram  of  the  fatigue  rig  Is  shown  In  Figure  3.1.  The 
test  rig  operates  on  the  four-square,  or  recirculating  torque,  principal.  As 
seen  In  Figure  3.1,  the  slave  gear  provides  the  recirculating  torque  by  virtue 
of  the  hydraulic  pressure  acting  on  the  load  vanes.  The  tooth  load  Is  con¬ 
trolled  by  varying  the  amount  of  hydraulic  pressure  applied.  Because  of  this 
recirculating  torque  principal,  the  drive  motor  provides  only  the  power 
required  to  overcome  the  frictional  losses  at  the  desired  running  speed.  The 
fatigue  rig  Is  capable  of  providing  75  kW  (100  hp)  at  the  designed  operating 
speed  of  10,000  rpm. 

3.1  TEST  GEARS 

The  gears  tested  using  the  fatigue  test  apparatus  were  all  the  same  type 
and  were  subjected  to  Identical  loading  conditions.  All  of  the  gears  were  made 
of  AISI  9310  steel  and  were  manufactured  to  AGMA  class  13,  aircraft  quality, 
tolerances.  The  test  gears  have  28  teeth,  a  pitch  diameter  of  88.9  mm 
(3.50  In.),  a  pressure  angle  of  20  degrees,  and  a  tooth  face  width  of  6.35  mm 
(C25  In.).  The  gears  were  loaded  to  74.6  Nm  (660  In  lbs),  which  resulted  In 
a  pltcMlne  maximum  hertz  stress  of  1.71  GPa  (248  ks  1  > ,  at  an  operating  speed 
of  10,000  rpm.  This  represents  a  load  of  almost  two  times  the  normal  design 
load  for  these  gears.  This  testing  procedure  was  used  In  order  to  obtain  fail¬ 
ures  within  a  reasonable  amount  of  test  time.  Table  I  lists  the  basic  run 
parameters  along  with  the  number  of  hours  to  failure  and  the  overall  mode  of 
failure  for  each  run.  As  seen  In  this  table,  the  only  run  parameters  that  dif¬ 
fer  between  the  runs  are  the  surface  treatment,  and  the  type  of  lubricant. 

Three  different  surface  treatments  were  encountered  In  the  eleven  runs 
recorded.  The  gears  In  runs  1  through  7  all  used  the  standard  surface  treat¬ 
ment,  which  Is  a  final  grinding  operation  to  a  surface  finish  of  356  mm  rms 
(14  in.  ,r. ) .  The  gears  In  runs  8  through  11  were  shot  peened  and  then  honed. 
The  main  ason  for  shot  peenlng  Is  to  create  subsurface  residual  stresses  that 
Improve  :ne  pitting  fatigue  life  of  a  gear.  The  gears  In  runs  8  through  10 
used  the  SPH  method,  a  high  Intensity  shot  peenlng  method  designed  to  produce 
high  subsurface  residual  stresses.  The  gears  In  run  11  used  the  SPL  method,  a 
low  Intensity  shot  peenlng  method  designed  to  produce  low  subsurface  residual 
stresses.  The  basic  differences  In  life  and  failure  modes,  as  seen  In  Table  I, 
are  primarily  due  to  the  different  lubricants  used. 
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TABLE  I:  Test  Run  Descriptions 


RUN  ft 

GEAR  SURFACE 
TREATMENT 

LUBRICATION 

TOTAL  LIFE 
(Hours) 

BASIC 

FAILURE  MODE 

1 

Standard 

Type 

A  Oil 

49 

Heavy  Wear 
and  Scoring 

2 

Standard 

Type 

A  Oil 

36 

Heavy  Wear 
and  Scoring 

3 

Standard 

Type 

A  Oil 

9 

Broke  Tooth 

4 

Standard 

Type 

A  Oil 

79 

Heavy  Wear 
and  Scoring 

5 

Standard 

Type 

A  Oil 

39 

Heavy  Wear 
and  Scoring 

6 

Standard 

Type 

A  Oil 

54 

Heavy  Wear 
and  Scoring 

7 

Standard 

Type 

A  Oil 

50 

Broke  Tooth 

8 

SPH 

Type 

B  Oil 

520 

No  Failure 
(Single  Pits) 

9 

SPH 

Type 

B  Oil 

245 

Single  Pits 

10 

SPH 

Type 

B  Oil 

339 

Distributed 

Pits 

11 

SPL 

Type 

B  Oil 

100 

Distributed 

Pits 
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The  two  different  oils  used  were  classified  type  A  oil  and  type  B  oil. 

The  type  A  oil  did  not  appear  to  have  sufficient  additives  to  provide  a  good 
Elasto-Hydro  Dynamic  (EHD)  film  thickness  between  the  gear  surfaces  under 
extreme  pressure.  Without  an  adequate  EHD  film  thickness  the  tooth  surfaces 
obtain  metal  to  metal  contact,  causing  severe  surface  wear  In  a  relatively 
short  period  of  time.  This  Is  especially  evident  by  the  fact  that  the  runs 
with  the  type  A  oil,  runs  1  through  7,  experienced  heavy  wear  and  surface  scor¬ 
ing,  and  an  average  life  of  only  63  hours.  Comparatively,  runs  8  through  11, 
which  used  type  B  oil,  experienced  an  average  life  of  300  hours,  and  subsur¬ 
face  fatigue  failures  In  the  form  of  pitting.  The  type  B  oil  was  a  synthetic 
paraffinic  oil  with  five  volume  percent  of  an  extreme-pressure  (EP)  additive. 
This  EP  additive  contains  sulfur  and  phosphorus  to  enable  the  oil  to  keep  an 
adequate  EHD  film  thickness,  minimizing  metal  to  metal  contact,  even  under 
extreme  pressure. 

The  failure  modes  experienced  by  the  various  gears  can  be  categorized  Into 
four  basic  damage  related  groups.  These  groups  are:  1)  Heavy  wear  and  scor¬ 
ing,  2)  Tooth  breakage,  3)  Single  pits,  and  4)  Distributed  pitting.  The  groups 
are  based  on  the  types  and  magnitude  of  damage  found  on  the  gears  at  the  end  of 
their  runs.  As  Indicated  In  Table  I,  run  8  was  not  classified  as  failed 
because  It  ran  the  maximum  of  520  hours  without  exceeding  the  rig's  vibration 
limit.  The  point  of  failure  on  the  fatigue  rig  was  determined  by  an  overall 
vibration  level  from  an  accelerometer  mounted  on  the  rig.  When  this  level 
reached  the  preset  threshold,  the  rig  shutdown  and  the  run  was  classified  as 
failed.  Because  the  gears  of  run  8  experienced  damage  similar  to  another  run. 
It  is  included  In  that  run's  group.  A  more  detailed  description  Is  given  In 
the  following  sections. 

HEAVY  WEAR  AND  SCORING: 

Runs  1,  2,  and  4  through  7  all  experienced  heavy  wear  and  scoring  uni¬ 
formly  over  all  the  gear  teeth.  Figures  3.2  and  3.3  Illustrate  the  wear  pat¬ 
terns  observed.  As  seen  In  these  figures,  the  driven  gear  teeth  experienced 
heavy  wear  and  scoring  at  the  tip  portion  of  the  tooth,  a  band  of  light  pitting 
at  the  double  to  single  tooth  contact  point,  and  heavy  wear  In  the  root  region 
of  the  tooth.  The  driver  gear  teeth  experienced  similar  wear  patterns  to  the 
driven  gear  teeth,  with  the  exception  that  the  top  portion  of  the  tooth  exper¬ 
ienced  lighter  wear  and  scoring.  Run  3  experienced  the  same  wear  patterns, 
although  to  a  lesser  degree  due  to  the  premature  tooth  failure  at  9  hours  Into 
the  test.  The  tooth  surfaces  appeared  extremely  course,  especially  In  the 
regions  of  heavy  wear  and  scoring. 

TOOTH  BREAKAGE 

Runs  3  and  7  experienced  tooth  breakage  during  the  fatigue  test  conducted. 
One  tooth  broke  during  run  3,  causing  the  rig  to  shut  down.  A  photograph  of 
the  fracture  surface  on  the  driven  gear  In  run  3  Is  shown  In  figure  3.4.  No 
detailed  tooth  crack  propagation  analysis  was  done  on  this  tooth,  however  by 
visual  Inspection,  there  appeared  to  be  only  three  distinct  bands  In  the  frac¬ 
ture  area,  with  little  evidence  of  fracture  surface  smoothing.  This  Indicates 
that  the  time  between  crack  formation  and  tooth  fracture  was  relatively  short. 
Run  7  also  experienced  tooth  fracture;  however,  a  failure  In  the  automatic  shut 
down  mechanism  did  not  stop  the  rig  soon  enough  to  prevent  extensive  damage  to 
both  gears.  Thus  no  tooth  fracture  observations  could  be  performed  on  run  7 
gears. 
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SINGLE  PITS 

Runs  8  and  9  experienced  similar  damage  In  the  form  of  single  pits. 
Figures  3.5  and  3.6  Illustrate  the  wear  pattern  found  on  the  gears  In  run  8  at 
the  end  of  the  run.  As  seen  In  these  figures  all  but  one  of  the  teeth  on  the 
driver  and  driven  gear  appear  to  have  no  sign  of  wear  or  pitting.  A  moderate 
pit  appears  on  one  tooth  on  the  driver  gear,  and  a  larger  pH  appears  on  one 
tooth  of  the  driven  gear.  Figures  3.7  and  3.8  Illustrate  the  condition  of  the 
gears  In  run  9  at  the  point  of  failure.  As  seen  In  these  figures,  run  9 
appears  to  have  a  similar  failure  pattern  as  found  In  run  8,  with  one  excep¬ 
tion.  The  large  pit  on  the  driven  gear  In  run  9,  although  roughly  the  same 
size  as  the  large  pit  In  run  8,  has  It's  major  axis  orientated  across  the  gear 
tooth  rather  than  along  the  tooth  surface,  as  found  In  run  8.  The  defects 
found  on  the  gears  In  runs  8  and  9  are  characterized  as  single  tooth  defects. 

DISTRIBUTED  PITTING 

Runs  10  end  11  experienced  similar  damage  In  the  form  of  distributed  pit¬ 
ting.  Figures  3.9  and  3.10  Illustrate  the  wear  pattern  present  on  the  gear 
teeth  In  run  10.  As  seen  In  these  figures,  the  driver  gear  exhibits  a  wear 
pattern  of  pitting  bands  with  various  Intensities  over  60  percent  of  the  teeth. 
The  driven  gear  displayed  a  similar  pattern  over  a  smaller  region  of  teeth. 
Figures  3.11  and  3.12  Illustrate  the  wear  pattern  present  on  the  gear  teeth  of 
the  gears  used  In  run  11.  As  seen  In  these  figures,  the  distributed  pitting 
bands  appear  In  a  similar  fashion  to  those  In  run  10,  although  to  a  lesser 
extent. 

3.2  INSTRUMENTATION  SET-UP 

Figure  3.13  Illustrates  the  test  Instrumentation  set-up  used  In  this 
study.  As  seen  In  this  figure,  two  accelerometers,  PCB  Model  303A,  were 
mounted  on  the  rig  close  to  the  test  gear  mesh.  The  vibration  signal  from 
these  accelerometers  were  recorded  on  a  high  precision  tape  recorder,  Sangamo 
Weston  Sabre  VII,  along  with  a  once  per  revolution  signal  and  a  time  code  sig¬ 
nal.  A  one  minute  signature  was  collected  every  three  hours  by  using  a  timer 
to  control  recording  time.  The  once  per  revolution  signal  was  provided  by  a 
photon  sensor  that  produced  a  narrow  time  pulse  (.202  msec)  for  each  revolution 
of  the  shaft.  This  signal  was  used  for  time  synchronous  averaging,  and  for 
determining  the  actual  rotational  speed.  The  time  code  signal  was  used  during 
tape  playback  to  distinguish  between  the  separate  data  points,  and  to  provide 
the  exact  time  and  day  for  each  data  point  recorded. 

After  the  runs  were  recorded,  the  data  was  then  analyzed  by  replaying  the 
tape  Into  a  single  channel  dynamic  signal  analyzer,  with  a  dynamic  range  of 
80  dB,  and  transferring  It  to  a  personal  computer.  The  raw  data  was  digitized 
and  averaged  by  the  analyzer  and  sent  via  a  general  purpose  Interface  bus  to  an 
IBM  compatible  personal  computer.  The  PC  used  a  standard  IEEE-488  Interface 
board  to  control  the  analyzer  and  download  the  digitized  data  to  DOS  files  on 
the  PC's  disc  drive.  The  data  transferred  was  the  averaged  amplitude  and  phase 
portions  of  the  Fourier  transform.  The  data  was  then  analyzed  by  using  several 
FORTRAN  programs,  each  based  on  the  various  prediction  techniques  investigated. 

The  primary  mesh  frequency  was  an  Important  factor  In  choosing  and  setting 
the  test  Instrumentation.  For  the  28  teeth  test  gears  operating  at  10,000  rpm, 
the  primary  and  first  harmonic  of  the  mesh  frequency  occurs  at  approximately 
4700  Hz  and  9400  Hz,  respectively.  The  accelerometers  were  chosen  with  a  fre¬ 
quency  range  of  0  to  10,000  Hz  to  capture  the  mesh  frequencies.  The  recorder 
tape  speed  was  set  at  762  mm/sec  (30  In/sec)  with  wideband  group  I  FM  recording 
electronics  to  give  a  recording  frequency  range  of  0  to  20  000  Hz,  with  a 
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dynamic  range  of  50  dB.  FM  recording  was  used  to  give  good  amplitude  sta¬ 
bility,  Independent  of  the  storage  characteristics  of  the  recording  tape. 

Using  a  tape  reel  with  a  2800  m  length  (9200  ft),  one  reel  translated  to 
184  hours  of  coverage,  or  7.7  days.  This  was  considered  a  reasonable  time 
Interval,  since  some  of  the  tests  ran  continuously  for  nearly  21  days.  The 
analyzer  was  set  up  to  time  average  the  signal  over  40  averages  for  each  data 
point,  using  the  shaft  synchronous  signal  as  the  Input  trigger.  The  FFT  was 
performed  over  the  range  0  to  12,500  Hz. 

During  the  first  stages  of  the  data  analysis  process,  It  was  discovered 
that  the  side  accelerometer  provided  the  same  Information  as  the  top  accelerom¬ 
eter.  In  fact,  the  meshing  frequency  components  In  the  side  accelerometer  sig¬ 
nal  were  not  as  dominant  as  In  the  top  accelerometer  signal.  This  Is  probably 
due  to  the  fact  that  the  top  accelerometer  Is  orientated  to  sense  the  vertical 
vibrations  of  the  box.  This  Is  coincident  with  the  direction  of  the  primary 
component  of  the  tooth  mesh  loads,  resulting  from  the  side  by  side  arrangement 
of  the  test  gears.  Thus,  the  vibration  data  from  the  top  accelerometer  was 
used  primarily  In  this  analysis,  with  the  side  accelerometer  data  used  as  a 
backup. 

3.3  FREQUENCY  RESPONSE  MEASUREMENTS  OF  FATIGUE  RIG 

During  the  Initial  stages  of  the  data  analysis,  It  was  discovered  that 
the  recorded  time  signals  contained  frequency  data  not  coherent  with  the  gear 
meshing  frequencies.  Figures  3.14  through  3.17  Illustrate  the  presence  of  the 

extraneous  frequencies  as  run  7  progressed  from  Initiation  to  failure.  This 

phenomenon  was  found  In  all  of  the  runs  recorded.  As  seen  In  these  figures, 
certain  frequencies,  not  coincident  with  the  gear  mesh  frequencies  or  asso¬ 
ciated  sidebands,  Increase  drastically  through  the  run  cycle.  The  primary  mesh 
frequency,  highlighted  by  an  arrow  In  the  figures,  Is  In  most  cases  lower  In 
magnitude  than  the  extraneous  frequencies.  The  amplitudes  near  3000  and 

6000  Hz  appear  to  be  the  most  dominant.  If  the  gear  failure  prediction  tech¬ 

niques  were  applied  to  the  signals  given  In  Figures  3.14  through  3.17,  results 
not  indicative  of  the  actual  gear  condition  would  dominate.  Therefore,  it  was 
necessary  to  Investigate  these  extraneous  frequency  components  and  determine 
their  origin.  If  these  components  could  be  attributed  to  a  known  source,  they 
then  could  be  removed  from  the  signal  with  confidence. 

A  number  of  possible  sources  of  the  extraneous  frequency  components  were 
investigated.  Because  the  signal  is  synchronous  time  averaged  with  the  gear 
shaft  rotation,  the  only  possible  periodic  sources  coincident  with  the  shaft 
speed  are  the  roller  bearings,  ball  bearings,  and  the  slave  gear.  The  location 
of  these  components  are  given  In  Figure  3.1.  Based  on  the  assumptions  of  a 
single  defect  and  that  the  rolling  elements  remain  in  rolling  contact,  bearing 
defect  frequencies  were  calculated  using  existing  equations.  ClU  These  fre¬ 
quencies  covered  outer  race,  inner  race,  and  rolling  element  defects,  and  fell 
within  the  64  to  1400  Hz  range.  These  explain  some  lower  frequency  components, 
but  they  were  too  low  for  the  dominant  frequencies  In  question.  The  slave  gear 
primary  mesh  frequency  was  found  to  be  6084  Hz,  which  is  near  some  of  the 
extraneous  frequencies,  but  it  was  not  exactly  coincident.  The  next  step  In 
Identifying  the  frequencies  In  question  was  to  determine  the  modal  parameters 
of  the  rig  using  frequency  response  measurements. 

Frequency  response  measurements  performed  on  the  fatigue  rig  Illustrated 
how  the  dynamics  of  a  transfer  path  can  critically  alter  a  signal.  The  fre¬ 
quency  response  measurements  were  taken  at  two  locations  on  the  box  with  the 
reference  point  being  the  top  accelerometer  location,  as  Illustrated  In 
Figure  3.18.  For  these  tests,  the  impact  method  was  used.  A  modal  hammer  with 
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a  steel  tip  and  force  transducer  was  used  as  the  Input  to  the  system,  with  the 
top  accelerometer  serving  as  the  output  signal.  The  force  transducer  and 
accelerometer  were  Input  Into  a  two  channel  dynamic  signal  analyzer  capable  of 
calculating  the  frequency  response  function.  Twenty  averages  were  taken  for 
each  measurement.  The  FFT  of  the  Input  Impulse  Is  given  in  Figure  3.19.  As 
seen  In  this  figure,  the  Impulse  adequately  excites  the  frequencies  In 
question. 

Results  from  the  frequency  response  measurements  directly  correlated  the 
unexplained  frequency  components  to  the  path  modes.  The  frequency  response 
measurement  of  the  path  from  the  top  of  the  bearing  cover  to  the  top  accelerom¬ 
eter  location  Is  given  In  Figure  3.20.  The  coherence  function  for  this  meas¬ 
urement,  shown  In  Figure  3.21,  Indicates  the  output  Is  almost  totally  a  func¬ 
tion  of  the  Input  up  to  10,000  Hz.  Figures  3.22  and  3.23  give  the  frequency 
response  measurement  and  Its*  coherence  for  the  path  from  the  gear  shaft  to  the 
top  accelerometer  location.  As  seen  In  Figures  3.20  and  3.22,  the  vibration 
path  from  the  gear  to  the  top  accelerometer  drastically  alters  the  signal  due 
to  the  natural  frequencies  Inherent  In  the  box.  As  the  gears  experience  wear, 
the  meshing  energy  redistributes  from  the  primary  meshing  frequencies  to  a  wide 
range  of  frequencies  In  the  frequency  spectrum.  Nhen  these  frequencies  become 
coincident  with  the  natural  frequencies  associated  with  the  transfer  path  the 
corresponding  amplitudes  Increase.  This  phenomenon  Is  clearly  Illustrated  In 
the  3000  Hz  region  and  In  the  region  within  the  5000  to  8000  Hz  band  In 
Figures  3.14  through  3.17.  The  frequencies  In  these  regions  are  coincident 
with  the  natural  frequencies  associated  with  the  transfer  path,  as  evident  In 
Figures  3.22  and  3.23. 

Most  of  the  frequencies  In  question  were  found  to  be  a  direct  result  of 
the  structural  dynamics  of  the  transfer  path  and  as  such  are  not  desired  as 
inputs  for  the  various  failure  prediction  techniques.  Therefore  these  frequen¬ 
cies  were  removed  from  the  processed  signals  by  filtering  techniques.  The 
functions  Illustrated  In  Figures  3.20  and  3.22  were  used  to  determine  the  fre¬ 
quencies  to  remove.  Careful  attention  was  placed  on  eliminating  only  those 
frequencies  known  to  be  dominant,  and  to  allow  frequencies  adjacent  to  the  mesh 
frequencies  to  be  passed.  Fortunately,  the  primary  mesh  frequency  region  had 
no  dominant  natural  frequencies  present;  however  some  of  the  upper  sidebands 
were  eliminated  with  the  filtering.  All  of  the  analysis  methods  were  applied 
to  the  filtered  signal.  It  Is  acknowledged  that  some  of  the  desired  signal  was 
eliminated  with  the  filtering;  however,  It  Is  assumed  that  the  errors  associated 
with  the  filtering  process  are  less  than  those  associated  with  leaving  the 
dynamic  effects  of  the  transfer  path  In  the  signal.  The  frequency  bands  most 
commonly  filtered  out  were  between  400  and  1500  Hz,  2400  and  3700  Hz,  and  5200 
and  8400  Hz. 


22 


DRIVER  GEAR 


Uniform  Damage 


Tooth  Specific  Damage 


Non* 


KEY 


m 

□m 

ms 

m 


HEAVY  SCORING 
UGHT  SCORING 
HEAVY  WEAR 
UGHT  PITTING  BAND 


DRIVEN  GEAR 


Uniform  Damage 


Tooth  Specific  Damage 


All  Dlmenelons 
In  mm 


None 


figure  3.2  :  Illustration  of  Tooth  Da*age  Found  on  Gears  in  Runs  1, 
2,  and  4  through  6. 
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Figure  3.5  :  Illustration  of  Tooth  Carnage  Found  on  Gears  in  Run  8. 
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b)  driven  gear  tooth 

Figure  ;  Photographs  of  tooth  damage  found 
on  gears  in  Run  8 
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Figure  3.7  :  Illustration  of  Tooth  Damage  Found  on  Gears  in  Run  9. 
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a)  driver  gear  tooth 


b)  driven  gear  tooth 

Figure  3.8  :  Photographs  of  tooth  damage  found 
on  gears  in  Run  9 
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Figure  3.9  :  Illustration  of  Tooth  Damage  Found  on  Gears  In  Run  10. 
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b)  Tooth  #  8,  driver  gear 
Figure  3.10  :  Photographs  of  tooth  damage  found 
on  gears  in  Run  10 
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d)  Tooth  #  10 t  driven  gear 

Figure  3.10  :  Concluded. 
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Figure  3.11:  Illustration  of  Tooth  Damage  Found  on  Gears  in  Run  11. 


34 


a)  Tooth  #  14,  driver  gear 


b)  Tooth  tt  16,  driver  goar 
Figure  3.12  :  Photographs  of  tooth  damage  found 
on  gears  in  Run  11 


c)  Tooth  #  8,  driven  gear 


d)  Tooth  #  4,  driven  gear 
Figure  3.12  :  Concluded. 
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Figure  3.13:  Illustration  of  Test  Set-up. 
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3.14:  Frequency  Spectrum  of  Run  7  at  Start  of  Run. 
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Figure  3. IS:  Frequency  Spectrum  of  Run  7  After  32  Hours  of  Running. 
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jure  3.16  :  Frequency  Spectrum  of  Run  7  After  44  Hours  of  Running. 
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Figure  3.17  :  Frequency  Spectrum  of  Run  7  After  50  Hours  of  Running. 
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Figure  3.20  :  Magnitude  of  Frequency  Response  Function,  Test  1, 

Between  the  Top  Accel ercaiater  Location  and  the  Top  of 
the  Bearing  Housing. 
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Figure  3,21  :  Coherence  Function  of  Test  1  Measurements,  Between  the 
Top  Accelerometer  Location  and  the  Top  of  the  Bearing 
Housing, 
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igure  3.22 


Magnitude  of  Frequency  Response  Function,  Test  2, 
Between  the  Top  Accelerometer  Location  and  the  Top  of 
the  Gear  Shaft. 
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Figure  3.23  :  Coherence  Function  of  Test  2  Measurements,  Between  the 
Top  Accelerometer  Location  and  the  Top  of  the  Gear 
Shaft. 


CHAPTER  IV 

APPLICATION  AND  RESULTS 


The  main  portion  of  this  study  was  the  application  of  the  various  gear 
fault  methods  to  experimental  data.  Chapter  4  briefly  explains  the  various 
computer  programs  developed  that  apply  each  of  the  techniques.  Also,  this 
chapter  presents  and  discusses  the  results  of  these  techniques. 

4.0  PROGRAMS  DEVELOPED 

To  apply  the  various  predictive  techniques  to  the  experimental  data,  sev¬ 
eral  computer  programs  were  developed.  The  programs  are  written  In  FORTRAN, 
and  listed  In  the  appendices.  For  those  methods  requiring  Fourier  analysis, 
the  standard  FFT  algorithm  developed  by  Cooley  and  Tukey  [123  Is  Included  In 
each  routine.  The  computer  programs  use  the  data  files  created  using  the 
dynamic  signal  analyzer.  The  results  of  the  programs  are  stored  In  data  files 
that  can  be  plotted  using  commercially  available  routines. 

The  programs  use  the  various  equations  ar.d  theories  presented  In  Chapter 
II.  Program  EFMO.FOR,  listed  In  Appendix  A,  calculates  the  FMO  parameter  for 
each  time  data  point  In  the  run.  Program  TFM4.F0R,  listed  In  Appendix  B, 
applies  the  FM4  techniques  to  the  data  by  calculating  the  normalized  kurtosls 
and  standard  deviation  values  of  the  difference  file  for  each  time  data  point. 
Program  HI  LB. FOR,  listed  In  Appendix  C,  uses  the  Hilbert  transform  technique 
to  calculate  the  Instantaneous  phase  variations  as  a  function  of  shaft  posi¬ 
tion  for  each  time  data  point.  Program  PARAM.FOR,  listed  In  Appendix  D,  calcu¬ 
lates  the  crest  factor,  sideband  level  factor,  and  energy  ratio  for  each  time 
data  point. 

The  computer  programs  were  verified  by  Inputting  known  functions.  As  an 
example  of  this,  the  kurtosls  and  standard  deviation  routines  In  program 
TFM4.F0R  were  verified  by  Inputting  a  sine  wave  and  square  wave,  and  comparing 
the  results  with  exact  solutions.  Table  II  Illustrates  the  results  of  this 
verification  test.  The  square  wave  did  not  provide  as  high  a  correlation  as  the 
sine  wave,  mainly  because  the  Input  square  wave  contained  some  Irregularities. 


TABLE  II:  Program  TFM4.FQR  routines  verification 


toutlne  -  Input 

Measured  value 

Actual  Value 

X  Difference 

Standard  Deviation 

-  Sine  wave 

1.3472 

1.3435 

0.28  X 

-  Square  Wave 

1.9180 

1.9000 

0.95  X 

Normalized  Kurtosls 

-  Sine  Wave 

1.4896 

1.5000 

0.69  X 

-  Square  Wave 

1 .0404 

1.0000 

4.04  X 

4.1  RESULTS 


FMO  METHOD 

The  FMO  method  did  very  well  at  detecting  a  majority  of  the  heavy  wear 
and  scoring  damage  experienced  by  runs  1,  2,  and  4  through  6,  as  It  was 
designed  to  do.  Figures  4.1  through  4.5  plot  the  three  parameters:  a)  FMO, 
b)  Peak-to-peak  level  (the  FMO  numerator),  and  c)  Sum  of  mesh  amplitudes  (the 
FMO  denominator),  for  these  runs.  As  seen  In  these  figures,  all  of  the  runs, 
except  run  1,  exhibit  an  Increase  in  the  FMO  parameter,  plot  a,  with  Increasing 
run  time.  These  same  runs  also  show  a  decrease  In  the  primary  meshing  frequency 
and  second  harmonic  amplitude  sum,  plot  c,  with  Increasing  run  time.  It  Is 
Interesting  to  note  that  In  Figure  4.3,  run  4,  the  FMO  plot  Indicates  heavy 
damage  at  49  hours  Into  the  test,  due  to  a  decrease  In  meshing  frequency  ampli¬ 
tudes,  whereas  the  peak-to-peak  level,  plot  b,  starts  Increasing  at  55  hours 
and  peaks  at  64  hours.  Run  1  did  not  provide  good  FMO  trends  even  though  the 
gear  experienced  similar  damage. 

FMO  Is  designed  to  detect  major  tooth  faults  such  as  the  breakage  exper¬ 
ienced  during  runs  3  and  7.  Because  no  data  was  collected  at  the  time  of  break¬ 
age,  or  Immediately  after  it  occurred,  FMO  was  unable  to  be  applied  to  these 
runs. 

FMO  was  applied  to  runs  8  and  9,  which  experienced  single  tooth  pits,  even 
though  FMO  Is  not  designed  to  detect  single  tooth  faults.  As  seen  In  Figure 
4.6,  run  8  shows  a  gradual  overall  Increase  In  FMO  values  with  run  time,  and  a 
gradual  overall  decrease  In  meshing  frequency  amplitudes.  This  trend  In  run  8 
could  be  attributed  to  Its  long  run  time  of  520  hours,  which  may  have  resulted 
In  uniform  wear  not  observable  by  visual  Inspection.  Run  9  showed  no  logical 
trend  with  FMO. 

FMO  was  also  applied  to  runs  10  and  11,  where  the  gears  experienced  dis¬ 
tributed  pitting  damage.  As  seen  in  Figures  4.7  and  4.8,  the  FMO  values  for 
both  runs  Increase  near  the  end  of  the  run  time.  The  meshing  frequency  ampli¬ 
tudes  show  an  overall  decrease  with  Increasing  run  time,  Indicating  the  exist¬ 
ence  of  heavy  wear.  Since  FMO  does  not  respond  to  specific  tooth  damage,  It 
Is  reasonable  to  assume  that  the  pitting  occurred  over  enough  teeth  to  act  as 
a  uniform  wear  phenomenon,  and  thus  was  capable  of  being  detected  by  FMO.  It 
is  Interesting  to  note  that  for  run  10,  Figure  4.7,  the  FMO  value  Increased 
sharply  at  the  220  hours  point.  At  this  same  point  In  time,  the  meshing  fre¬ 
quencies  amplitude  sum  decreased  sharply,  possibly  indicating  the  time  when 
most  of  the  major  distributed  pitting  happened.  The  peak-to-peak  level  at 
this  period  contained  no  major  fluctuation,  only  a  gradual  Increase  In  level. 

FM4  METHOD 

The  normalized  kurtosls  parameter  was  unable  to  detect  the  single  tooth 
faults  found  In  the  test  gears.  Runs  8  and  9  experienced  specific  faults  In 
the  form  of  one  large  and  one  small  pit.  As  seen  In  Figures  4.9  and  4.10, 
these  defects  had  limited  effect  on  the  kurtosls  parameter.  The  fatigue  cracks 
that  resulted  In  broken  teeth  In  runs  3  and  7  were  also  undetected  by  the 
kurtosls  parameter,  as  seen  In  Figure  4.11  for  run  7.  The  normalized  kurtosls 
values  for  the  two  data  points  of  run  3,  not  plotted,  were  2.98  and  3.04.  It 
should  be  noted  that  the  data  points  collected  at  the  end  of  runs  3  and  7  were 
taken  3  and  2  1/2  hours,  respectively,  prior  to  the  point  of  tooth  fracture. 

Only  run  11  had  limited  success  using  the  normalized  kurtosls  parameter.  As 
seen  In  Figure  4.12,  a  normalized  kurtosls  value  of  5.3  was  registered  at 
88  hours  of  operation;  however,  the  value  decreased  to  approximately  3.0  at 
96  hours.  One  explanation  for  this  trend  Is  that  only  one  large  pit  was 
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present  at  the  88  hour  mark.  As  more  pits  formed,  the  number  of  peaks 
Increased  causing  the  normalized  kurtosls  value  to  decrease. 

The  standard  deviation  parameter,  although  used  In  FM4  for  single  tooth 
failure  detection,  proved  to  be  a  good  Indicator  of  heavy  wear.  Runs  1,  2  and 
4  through  6  all  experienced  heavy  wear  and  scoring.  The  standard  deviation 
plots  for  these  runs,  given  In  Figures  4.13  through  4.17,  all  show  clear  trends 
of  Increasing  values  near  the  end  of  each  run  with  only  minor  fluctuations. 

The  standard  deviation  plot  of  run  1,  Figure  4.13,  Indicates  a  definite  trend 
as  compared  to  the  FMO  plot  of  run  1,  which  did  not  respond  to  the  heavy  wear. 

Plots  of  the  squared  difference  signal  of  FM4  provided  some  general  Infor¬ 
mation  on  the  wear  condition  of  the  tooth  surfaces.  The  plots  of  the  squared 
difference  signal  are  presented  as  a  function  of  gear  rotation.  Two  plots  are 
given  In  each  figure.  The  top  graph  (a)  represents  the  data  collected  at  the 
beginning  of  the  run  and  the  bottom  graph  (b)  represents  the  data  collected  at 
the  end  of  the  run.  The  runs  that  experienced  heavy  wear  and  scoring  show  def¬ 
inite  Increases  In  the  difference  signal  between  the  start  and  iid  plots 
throughout  the  gear  rotation.  An  example  of  this  Is  given  In  Figure  4.18.  The 
specific  pits  of  runs  8  and  9  were  not  seen  in  their  squared  difference  signal 
plots.  The  wear  pattern  on  the  driver  gear  In  run  10  Is  reflected  In  the 
squared  difference  signal  given  In  Figure  4.19.  Run  11  experienced  similar 
wear  as  run  10;  however,  the  wear  pattern  was  not  easily  detected  In  Its 
squared  difference  signal  plot. 

HILBERT  TRANSFORM  METHOD 

The  Hilbert  transform  method  was  primarily  developed  to  detect  fatigue 
cracks  using  the  phase  modulation  function  described  In  Section  2.1.  Runs  3 
and  7  are  the  only  runs  that  experienced  tooth  fracture  due  to  probable  fatigue 
cracks.  Figure  4.20  plots  the  phase  modulation  function  for  the  two  data  time 
Intervals  of  run  3.  Several  phase  shifts  can  be  seen  In  the  first  plot,  2 
hours  Into  the  run;  however,  they  are  not  reflected  In  the  second  plot  repre¬ 
senting  the  last  data  point.  Nothing  In  the  last  data  point  plot  Indicates 
the  presence  of  a  fatigue  crack,  l.e.  no  large  phase  lags  present.  Figure  4.21 
plots  the  phase  modulation  function  for  the  last  two  data  time  Intervals  of 
run  7.  The  only  possible  Indication  of  a  fatigue  crack  starting  Is  the  phase 
lag  at  approximately  half  way  Into  the  shaft  rotation,  near  the  180  degrees 
point.  The  phase  lag  starts  during  the  second  to  last  time  Interval,  47  hours 
Into  the  run,  and  grows  to  the  size  seen  In  the  last  Interval,  50  hours  Into 
the  run.  More  data  after  this  last  data  point  Is  required  to  claim  with  any 
certainty  that  this  phase  lag  represents  an  accual  fatigue  crack.  Again  It 
must  be  noted  that  the  last  data  points  for  runs  3  and  7  were  taken  3  and  2 
1/2  hours,  respectively,  before  tooth  fracture. 

CREST  FACTOR  AND  SIDEBAND  LEVEL  FACTOR 

Both  the  crest  factor  and  sideband  level  factor  are  designed  to  respond 
to  signals  with  Impulsive  vibration  sources,  specifically  tooth  breakage. 

Runs  3  and  7  were  the  only  run^  that  experienced  tooth  fracture.  The  crest 
factor  and  sideband  level  factor  for  run  7  are  plotted  In  Figure  4.22.  Run  3 
was  not  plotted,  since  It  had  only  two  time  Intervals  recorded  before  fracture 
occurred.  As  seen  In  Figure  4.22,  neither  tho  crest  factor  (plot  a)  or  the 
sideband  level  factor  (plot  b)  display  any  Indication  that  a  tooth  fracture 
was  going  to  occur.  These  parameters  may  be  sensitive  to  tooth  breakage  only 
after  It  has  happened.  Unfortunately,  no  data  was  collected  during  or  after 
the  fracture  occurred. 
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ENERGY  RATIO 

The  energy  ratio  Is  designed  to  be  a  robust  Indicator  of  heavy  wear. 

Runs  1,  2,  and  4  through  6  experienced  heavy  wear  and  scoring.  The  energy 
ratio  graphs  of  these  runs  are  given  In  Figures  4.23  through  4.27.  As  seen  In 
these  figures,  the  energy  ratio  does  not  provide  as  good  an  Indication  of  wear 
as  the  FMO  and  the  standard  deviation  methods.  The  energy  ratio  graphs  for 
runs  10  and  11  are  given  In  Figure  4.28  and  4.29,  respectively.  Of  these  two 
runs,  only  run  10  had  an  Increase  In  Its*  energy  ratio  parameter. 
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Figure  4.2  :  Plot  of  FMO  Parameters  for  Run  2. 

a)  FHO  Parameter. 

b)  Numerator  of  FMO. 

c)  Denominator  of  FMO. 
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Figure  4.4  :  Plot  of  FMO  Parameters  for  Run  5 

a)  FMO  Parameter. 

b)  Numerator  of  FMO. 

c)  Denominator  of  FMO. 
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Plot  of  FMO  Parameters  for  Run  6 

a)  FMO  Parameter. 

b)  Numerator  of  FMO. 

c)  Denominator  of  FMO. 
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Figure  4.6  :  Plot  of  FNO  Parameters  for  Run  8. 

a)  FHO  Parameter. 

b)  Numerator  of  FHO. 

c)  Denominator  of  FHO. 
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Figure  4.9  :  Plot  of  Normalized  Kurtosis  (FN4)  for  Run  8. 
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Figure  4.10  :  Plot  of  Normalized  Kurtosis  (FN4)  for  Run  9. 
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Plot  of  Normal Ized  Kurtosls  (FM4)  for  Run  7. 
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Plot  of  Normalized  Kurtosls  (FM4)  for  Run  H 
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Figure  4.13  :  Plot  of  Standard  Deviation  of  Difference  Signal  (FM4) 
for  Run  1. 


Figure  4.14  :  Plot  of  Standard  Deviation  of  Difference  Signal  (FN4) 
for  Run  2. 
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Figure  4.17  :  Plot  of  Standard  Deviation  of  Difference  Signal  (FM4) 
for  Run  6. 
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Figure  4.18  :  Plot  of  Squared  Difference  Signal  For  Run  4. 

a)  At  Start  of  Run. 

b)  At  End  of  Run. 
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figure  4.19  :  Plot  of  Squared  Difference  Signal  for  Run  10. 

a)  At  Start  of  Run. 

b)  At  End  of  Run. 
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Figure  4.21 
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Plot  of  Phase  Modulation  Function  for  Run  7. 

a)  At  a  Run  Time  of  47  Hours. 

b)  At  a  Run  Time  of  SO  Hours. 
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Figure  4,22  :  Plot  of  Crest  Factor  and  Sideband  Level  for  Run  7. 

a)  Crest  Factor. 

b)  Sideband  Level  Factor. 
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Figure  4.23  : 


Energy  Ratio  Plot  for  Run  1. 
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Figure.  4.25  :  Energy  Ratio  Plot  for  Run  4. 


Figure  4.26  :  Energy  Ratio  Plot  for  Run  5. 
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CHAPTER  V 

CONCLUSIONS  AND  RECOMMENDATIONS 


The  type  and  extent  of  the  damage  found  on  the  eleven  test  gears  can  be 
classified  Into  four  major  failure  modes.  The  first  mode  can  be  described  as 
heavy  wear  and  scoring  (runs  1,2  and  4  through  6).  The  second  mode  Is  tooth 
breakage  (runs  3  and  7).  The  third  failure  mode  Is  a  result  of  single  pits 
(runs  8  and  9).  And  the  last  mode  Is  described  as  distributed  pitting  (runs 
10  and  11).  These  classifications  will  be  used  to  evaluate  the  overall  per¬ 
formance  of  the  various  methods. 

HEAVY  WEAR  AND  SCORING 

The  FMO  parameter  and  the  standard  deviation  of  the  difference  signal 
from  FM4  did  very  well  at  detecting  this  heavy  wear  condition.  Results  from 
both  techniques  support  the  theory  that  as  a  gear  wears  the  meshing  energy 
redistributes  from  the  meshing  frequencies  to  Its  sidebands  and  beyond. 

The  energy  ratio  method  does  not  predict  wear  as  well  as  the  FMO  and  the 
standard  deviation  methods.  The  most  probable  reason  for  this  Is  In  the  denom 
Inator  of  the  energy  ratio.  The  denominator  Is  the  standard  deviation  of  the 
regular  signal,  or  the  meshing  frequencies  and  their  first  order  sidebands. 

The  first  order  sidebands  should  not  be  considered  part  of  the  regular  signal 
for  this  parameter.  This  Is  because  they  Increase  In  amplitude  similar  to  the 
higher  order  sidebands  as  the  gear  wears. 

An  enhanced  method  using  a  combination  of  these  techniques  could  prove  to 
be  a  reliable  uniform  gear  wear  detection  technique.  One  such  method  could 
use  the  standard  deviation  of  the  difference  signal  divided  by  the  sum  of  the 
amplitudes  of  the  meshing  frequencies  (primary  and  first  harmonic).  Because 
wear  becomes  detectable  In  the  overall  vibration  level  only  when  it  reaches  a 
critical  stage,  a  robust  wear  indication  parameter  would  be  highly  useful. 

TOOTH  BREAKAGE 

The  fatigue  cracks  that  resulted  In  the  broken  teeth  in  runs  3  and  7  were 
no1^  detected  by  any  of  the  methods.  The  normalized  kurtosls  parameter  of  FM4, 
the  phase  modulation  function  from  the  Hilbert  transform  technique,  the  crest 
factor,  and  the  sideband  level  factor  are  all  conditioned  to  react  to  tooth 
cracks  to  varying  degrees;  however,  none  detected  any  fault  prior  to  tooth 
fracture.  The  most  probable  reason  for  this  can  be  due  to  the  high  speed  and 
high  loading  conditions  the  test  gears  are  subjected  to.  With  these  operating 
conditions,  the  time  elapsed  between  Initiation  of  the  fatigue  crack  and  even¬ 
tual  tooth  fracture  Is  probably  orders  of  magnitude  lower  than  the  three  hour 
data  Interval  time.  Since  the  last  time  records  collected  from  runs  3  and  7 
were  3  and  2  1/2  hours  before  actual  tooth  fracture,  respectively,  data 
necessary  to  indicate  the  fatigue  cracks  were  missed. 

SINGLE  PITS  (TOOTH  SPECIFIC  FAULTS) 

The  single  pits  were  not  detected  by  the  normalized  kurtosls  In  the  FM4 
technique.  The  normalized  kurtosls  of  the  difference  signal  Is  designed  to 
detect  tooth  specific  faults;  however,  no  Indication  of  these  defects  were 
recorded  with  this  parameter.  Although  they  physically  appear  large,  they  may 
have  not  affected  the  signal  to  the  point  detectable  by  FM4. 

DISTRIBUTED  PITTING 

The  distributed  pitting  damage  was  detected  with  the  the  FMO  parameter. 

It  Is  theorized  that  the  pitting  happened  over  enough  of  the  teeth  to  act  as  a 
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uniform  wear  phenomenon,  thus  becoming  detectable  to  FMQ.  The  square  of  the 
difference  signal  from  the  FM4  method  did  appear  to  reflect  the  actual  wear 
pattern  on  some  of  the  runs. 

RECOMMENDATIONS 

It  Is  recommended  that  any  future  tests  of  this  nature  be  Implemented  on 
an  on-line  basis  for  three  major  reasons.  The  first  stems  from  the  fact  that 
a  passive  system  cannot  call  attention  to  events  which  happen  faster  than  the 
time  between  the  data  time  Intervals.  A  continuous,  or  on-line,  system  can 
monitor  the  test  and  record  data  as  fast  as  the  electronics  will  allow  when 
certain  parameters  signal  a  significant  event  In  progress.  This  may  be  the 
only  way  to  provide  data  on  the  fatigue  cracks  that  lead  to  tooth  fracture  on 
this  rig.  Another  major  reason  Is  that  In  an  on-line  system  It  Is  possible  to 
have  the  system  stop  the  rig  when  a  parameter  detects  a  fault.  The  fault  can 
be  analyzed  to  determine  Its  extent  and  nature  at  the  exact  time  It  Is 
detected.  This  would  provide  an  excellent  means  of  verifying  and  linking  cer¬ 
tain  parameters  with  certain  faults  as  they  naturally  happen.  The  final  advan¬ 
tage  of  an  on-line  system  Is  the  ability  of  storing  data  In  digital  form, 
eliminating  the  FM  recorder.  In  this  test,  because  the  dynamic  range  of  the 
recorder  was  only  50  dB,  some  of  the  lower  energy  components  of  the  signal  were 
lost.  Using  the  analyzer  to  digitize  and  store  the  data  directly  would 
Increase  the  dynamic  range  to  80  dB,  thereby  capturing  more  of  the  lower  energy 
signals. 

It  Is  also  recommended  that  frequency  response  measurements  be  performed 
between  the  gear  shaft  and  the  proposed  measurement  sites.  The  dynamics  of  the 
transfer  path  In  this  study  were  found  to  significantly  affect  the  output  sig¬ 
nal,  especially  as  the  gear  began  to  wear.  An  understanding  of  the  structural 
dynamics  associated  with  the  transfer  path  Is  required  to  apply  prediction 
techniques  only  to  the  portion  of  the  signal  reflecting  the  actual  gear  mesh 
dynamics. 
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APPENDIX  A  -  SOURCE  CODE  FOR  EFMO.FOR 

OOOl  C*********************************^*******^********^**************C 


0002  C  C 
0003  C  C 
0004  C  PROGRAM:  EFMO.FOR  CREATED:  05-05-89  REVISED:  XX-XX-XXXX  C 
0005  C  C 
0006  C  PROGRAMMER:  J.J.  ZAKRAJSEK  C 
0007  C  C 
0008  C  THIS  PROGRAM  COMPUTES  AN  ENHANCED  VERSION  OF  THE  COURSE  FAULT  C 
0009  C  DETECTION  PARAMETER  (EFMO) .  IT  IS  THE  BASIC  FMO  PARAMETER,  C 
0010  C  AS  PROPOSED  BY  STEWART  IN  HIS  1977  PAPER,  HOWEVER,  THE  PEAK  C 
0011  C  TO  PEAK  VALUE  USED  IN  EFMO  IS  THE  PEAK  TO  PEAK  OF  THE  C 
0012  C  RECONSTRUCTED  TIME  SIGNAL  AFTER  BAND  WIDTH  FILTERING  IS  C 
0013  C  PERFORMED.  THE  RECONSTRUCTED  TIME  SIGNAL  USED  IS  THE  "FTIME"  C 
0014  C  FILE  CREATED  AFTER  A-B-C  FILTERING  DURING  FM4  CALCULATIONS  IN  C 
0015  C  PROGRAM  TFM4.FOR.  THE  AMPLITUDES  AT  THE  MESH  FREQUENCY  AND  ITS  C 
0016  C  HARMONIC  ARE  RETRIEVED  FROM  PREVIOUS  FMO  CALCULATIONS.  C 
0017  C  C 
0018  C  C 


0019  C****************************************-*******'******'*****'>|f*,***'*****C 

0020  C 

0021  DIMENSION  TIME(50) , ECFDP( 50) , PKTPK(50) , ASUM(50) , FSIG(220) 

0022  CHARACTERS  FMOF, EFMOF, LI ST2 .FTIMEF 

0023  C 

0024  C  ENTER  INITIAL  DATA  AND  FILE  NAMES 
0025  C 

0026  PRINT  *,  'ENTER  NUMBER  OF  DATA  POINTS' 

0027  READ(*,*)  K 

0028  PRINT  *,  'ENTER  OUTPUT  DATA  FILE  NAME  (\TXXE*FM0. DAT) ' 

0029  PRINT  *,  •f  *  -  FILTER  CONFIG  USED  IN  RECONSTRUCTED  FILTER)' 

0030  READ  <*  (A22) ' )  EFMOF 

0031  OPEN(7 , FILE-EFMOF, IOSTAT-IQERR) 

0032  IF  (IOERR  .GT.  0)  GOTO  999 

0033  PRINT  *,  'ENTER  FILE  NAMES  ARRAY  FILE  (\TXXLIST2 . FIL)  * 

0034  READ(*, ' (A22) ')  LIST2 

0035  OPEN ( 8 , FIT£-LIST2 , IOSTAT-IOERR) 

0036  I F( IOERR  .GT.  0)  GOTO  999 

0037  PRINT  *,  'ENTER  PREVIOUS  TFHO  FILE  (\TXXFMO.DAT)' 

0038  READ(*, ' (A22) ' )  FHOF 

0039  OPEN(9 , FI LE-FHOF, IOSTAT-IOERR) 

0040  IF( IOERR  .GT.  0)  GOTO  999 

0041  PRINT  *,  'ENTER  NUMBER  OF  TIME  DATA  POINTS  PER  CYCLE  (INTEGER)' 

0042  READ(*,*)  NCYC 

0043  C 

0044  C  START  TIHE  LOOPING,  PROCESSING  DATA  AT  EACH  TIME  INTERVAL 
0045  C 

0046  DO  300,  J-l.K 

0047  PRINT  *,  'ACCESSING  AND  PROCESSING  DATA  FILE  NO: ' ,J 

0048  READ(8 ,800)  FTIMEF 

0049  800  F0RMAT(AI8) 

0050  READ(9 ,900)  TIHK(J) ,ASUN(J) 

0051  900  FORMAT(FB . 2 , 44X , F12 .8) 

0052  0PEN(4 , FILE-FTIMEF, IOSTAT-IOERR) 

0053  IF  (IOERR  .GT.  0)  GOTO  999 

0054  C 

0055  C  READ  IN  AND  FIND  THE  MAX  AND  MIN  VALUES  OF  THE  FILTERED  TIME  RECORD 
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0056  C 

0057  AMAX  -  0.0 

0058  BMAX  -  0.0 

0059  DO  100,  I-l.NCYC 

0060  READ(4,400)  FSIG(I) 

0061  400  FORMAT ( 20X , F16 . 8 ) 

0062  IF(  FSIG(I)  .GT.  AMAX  )  THEN 

0063  AMAX  -  FSIG(I) 

0064  END IF 

0065  IF(  FSIG(I)  .LT.  BMAX  )  THEN 

0066  BMAX  -  FSIG(I) 

0067  ENDIF 

0068  100  CONTINUE 

0069  C 

0070  C  DETERMINE  PK  TO  PK  OF  FILTERED  SIGNAL,  AND  ENHANCED  FMO  VALUE 
0071  C 

0072  PKTPK(J)  -  AMAX  -  BMAX 

0073  ECFDP(J)  -  PKTPK ( J ) /ASUM ( J ) 

0074  C 

0075  C  CLOSE  LOOP  DEPENDANT  FILE,  AND  WRITE  CALAULTED  DATA 
0076  C 

0077  CLOSE (4 , IOSTAT-IOERR) 

0078  IF  (IOERR  .GT.  0)  GOTO  999 

0079  WRITE  (7,700)  TIME(J) ,ECFDP(J) ,PKTPK(J) 

0080  700  FORMAT ( F8 . 2 , F10 . 4 , F10 . 4) 

0081  300  CONTINUE 

0082  C 

0083  PRINT  *,  'EFMO  VALUES  STORED  IN  FILE ; * , EFMOF 

0084  STOP 

0085  999  PRINT*,  '  ERROR  NO . \IOERR 

0086  END 
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APPENDIX  B  -  SOURCE  CODE  FOR  TFM4.F0R 

C********************************************************************C 


C  C 
C  C 
C  PROGRAM:  TFM4 . FOR  REVISION:  NEW  C 
C  REVISION  DATE:  NEW  C 
C  PROGRAMMER:  J.J.  Zakrajsek  DATE  CREATED:  04-13-89  C 
C  C 
C  THIS  PROGRAM  COMPUTES  THE  FAULT  DETECTION  PARAMETERS  AND  THE  C 
C  RELATED  MATHEMATICS  FOR  THE  FM4  DIAGNOSTIC  THEORY  PROPOSED  C 
C  IN  STEWART'S  1977  PAPER  "  SOME  USEFULL  DATA  ANALYSIS  C 
C  TECHNIQUES  FOR  GEARBOX  DIAGNOSTICS . •  THESE  INCLUDE  PARAMETERS  C 
C  SUCH  AS  STANDARD  DEVIATION  AND  KURTOSIS.  C 
C  C 
C  C 


C*^^*************^*,iWf****^**^**************ilWr**'******'***********'*C 

C 

DIMENSION  TIME(50) ,DYO(220) ,DSYO(220) ,RMAG(401) 

DIMENSION  FREQ(401) ,FMAGRMS(401),FMAG(401),PHAS(401),FFMAG(401) 
DIMENSION  TO(1024) ,YO(1024) , FY0(1024) ,RFYO(1024) 

CHARACTER*22  TIMEF , FM4F, LIST1 , LIST2 , AMPLF , PHASEF . ATIME 
CHARACTER*! 8  FTIME , RTIME , NTIME, STIME 
COMPLEX  X<1024) 

C 

C  ENTER  GENERAL  PARAMETERS  AND  FILE  NAMES 
C 

PRINT  *,  'ENTER  NUMBER  OF  DATA  POINTS  (INTEGER)' 

READ<*,*)  K 

PRINT  *,  'ENTER  MESHING  FREQUENCY.  (XXXX.XX)' 

READ(*,*)  PRIM 

PRINT  *,  'ENTER  ROTATIONAL  SPEED  IN  REV/SEC  (Hz)' 

READ(*,*)  SPD 

PRINT  *,  'ENTER  NUMBER  OF  TIME  DATA  POINTS  PER  CYCLE' 

PRINT  *,  '(  INTEGER,  220  MAXIMUM  )' 

READ  (*,*)  NCYC 
C 

C  DEFINITION  OF  OPTIONAL  FILTER  BANDS  A,  B,  AND  C 
C 

PRINT  *,  '3  BAND  WIDTH  DIGITAL  FILTERS  CAN  BE  SELECTED.' 

PRINT  *,  ’DATA  EXISTING  WITHIN  THE  THREE  BANDS  WILL  BE  FILTERED' 
PRINT  *.  'OUT  OF  THE  ANALYSIS  BEFORE  IT  BEGINS  FM4  CALCULATIONS' 
PRINT  *.  ' 

PRINT  *,  'TO  USE  FILTERS  ENTER  (1)’ 

PRINT  *,  'TO  SKIP  FILTERS  ENTER  (0)' 

C 

READ  (*,+)  NN 
IF  (NN  .LT.  I)  GOTO  123 
C 

PRINT  *.  'ENTER  STARTING  FREQUENCY  OF  FILTER  A' 

READ(*,*)  AFILTS 

PRINT  *,  '***************  ENTER  END  FREQUENCY  OF  FILTER  A' 
READ(*,*)  AFILTE 

PRINT  *,  'ENTER  STARTING  FREQUENCY  OF  FILTER  B' 

READ(*,*)  BFILTS 

PRINT  *,  ' ***************  ENTER  END  FREQUENCY  OF  FILTER  B' 
READ(*,*)  BFILTE 
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PRINT  *,  ' ENTER  STARTING  FREQUENCY  OF  FILTER  C' 

READ(*,*)  CFILTS 

PRINT  *,  ' ***************  ENTER  END  FREQUENCY  OF  FILTER  C' 
READ(*,*)  CFILTE 
123  CONTINUE 
C 

C  INPUT  PRIMARY  FILE  NAMES  FILES 
C 

PRINT  *,  'ENTER  TXXTIME. FIL  INPUT  FILE  NAME' 

READ(*, ' (A22) ' )  TIMEF 

PRINT  *,  'ENTER  TXXFM4*.FIL  OUTPUT  FILE  NAME' 

PRINT  *,  '(  *  -  FILTERS  A-B-C  VERSION)' 

READ(*, ' (A22) ' )  FM4F 

PRINT  *,  'ENTER  TXXLIST1.FIL  INPUT  FILES  FILENAME' 
READ(*,'(A22)')  LIST1 

PRINT  *,  'ENTER  TXXLIST2 . FIL  OUTPUT  FILES  FILENAME* 
READ(*, ' (A22) ' )  LIST2 
C 

C  OPEN  PRIMARY  FILES 
C 

OPEN ( 1 , FI LE-TIMEF, IOSTAT-IOERR) 

IFUOERR  .GT.  0)  GOTO  999 
OPEN< 2 .FILE-FM4F, IOSTAT-IOERR) 

IF<IOERR  .GT.  0)  GOTO  999 
OPEN ( 3 , FILE-LIST1 , IOSTAT-IOERR) 

IFUOERR  .GT.  0)  GOTO  999 
0FEN(4 , FILE-LIST2 , IOSTAT-IOERR) 

IFUOERR  .GT.  0)  GOTO  999 
C 

C  STARTING  POINT  OF  FM4  ANALYSIS  LOOP  FOR  EACH  DATA  POINT 
C 

DO  111,  J-i.K 

PRINT  *.  'PROCESSING  DATA  POINT  NO.  : » , J 
READ( 1,100)  TIME(J) 

100  FORMAT (F10. 2) 

READ(3,300)  AMPLF, PHASEF.ATIHE 
300  PORHAT(A22,A22,A22) 

READ (4, 400)  FTIHE , RTIME , NTIHE , STIME 
400  F0RMAT(A18,A18,A18,A18) 

C 

C  OPEN  LOOP  DEPENDANT  FILES 
C 

OPEN ( 5 , FI LE-AHPLF , IOSTAT-IOERR) 

IFUOERR  .GT.  0)  GOTO  999 
OPEN ( 6 , FI LE-PHASEF , IOSTAT-IOERR) 

IFUOERR  .GT.  0)  GOTO  999 
OPEN<  7 , FI LE-ATIHE, IOSTAT-IOERR) 

IFUOERR  .GT.  0)  GOTO  999 
OPEN (8 , FI LE-FTIHE, IOSTAT-IOERR) 

IFUOERR  .GT.  0)  GOTO  999 
OPEN (9 , F1LE-RTIHE, IOSTAT-IOERR) 

IFUOERR  .GT.  0)  GOTO  999 
OPEN(iO,FILE-NTIHE, IOSTAT-IOERR) 

IFUOERR  .GT.  0)  GOTO  999 
OPEN(il , FILE-STINE, IOSTAT-IOERR) 
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0111  IF(I0ERR  .GT.  0)  GOTO  999 

0112  C 

0113  C  READ  IN  PHASE  AND  FREQUENCY  FILES 
0114  C 

0115  DO  222,  1-1,401 

0116  READ(5,500)  FREQ( I ) , FMAGRMS ( I ) 

0117  500  FORMAT(4X,F12.3,20X,F16.8) 

0118  FHAG(I)  -  FMAGRMS ( I )  *  (2. **.5) 

0119  READ(6,600)  PHAS(I) 

0120  600  FORMAT (2 OX, Fl 2 . 3) 

0121  222  CONTINUE 

0122  C 

0123  C  PERFORM  IFFT  ON  ORIGINAL  SIGNAL 
0124  C 

0125  CALL  IFOUR(FREQ , FMAG , PHAS , TO, YO) 

0126  C 

0127  C  PRINT  TO  FILE,  2  CYCLES  OF  ORIGINAL  TIME  SIGNAL. 

0128  C 

0129  NCYC2  -  NCYC  +  NCYC 

0130  DO  333,  I-1.NCYC2 

0131  WRITE(7 , 700)  TO(I),YO(I) 

0132  700  FORMAT(F16 . 8 ,4X,F16 .8) 

0133  333  CONTINUE 

0134  C 

0135  C  PERFORM  FILTERING  USING  FILTERS  A,  B,  C 
0136  C 

0137  IF<NN  ,LT.  1)  GOTO  234 

0138  CALL  FILTER( FREQ , FMAG , FFMAC , AFILTS , AFILTE , BFILTS , BFILTE , 

0139  *CFILTS , CFILTE) 

0140  C 

0141  C  PERFORM  IFFT  ON  FILTERED  SIGNAL 
0142  C 

0143  CALL  I FOUR( FREQ , FFMAG , PHAS , TO , FYO) 

0144  234  CONTINUE 

0145  IF(NN  ,LT.  1)  THEN 

0146  DO  345,  1-1,401 

0147  FFKAG(I)  -  FMAG(I) 

0148  345  CONTINUE 

0149  DO  456,  I-1.NCYC2 

0150  FYO(I)  -  YOU) 

0151  456  CONTINUE 

0152  ENDIF 

0153  C 

0154  C  PERFORM  FILTERING  TO  FK4  SPECIFICATIONS 
0155  C 

0156  CALL  FFM4 (FREQ, FFMAG, PRIH.SPD.RMAG) 

0157  C 

0158  C  PERFORM  IFFT  ON  FN4  FILTERED  SIGNAL 
0159  C  (USES  ABC  FILTERED  SIGNAL  AS  BASIS) 

0160  C 

0161  GALL  IFOUR<  FREQ , RMAC , PHAS ,T0 , RFYO) 

0162  C 

0163  C  DETERMINE  DIFFERENCE  SIGNAL  <  ORIGINAL  •  REGULAR  COMPONENTS) 
0164  C 

0165  DO  444,  I-i,NCYC 
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DYO(I)  -  FYO(I)  -  RFYO(I) 

444  CONTINUE 
C 

C  DETERMINE  SQUARED  DIFFERENCE  SIGNAL 
C 

DO  555,  I-l.NCYC 
DSYO(I)  -  DY0(I)**2. 

555  CONTINUE 
C 

C  DETERMINE  AVERAGE  VALUE  OF  DIFFERENCE  SIGNAL 

DSUM  -  0.0 
DO  666,  I-l.NCYC 
DSUM  *■  DSUM  +  DYO(I) 

666  CONTINUE 

DBAR  -  DSUM/FLOAT(NCYC) 

C 

C  DETERMINE  STANDARD  DEVIATION  OF  NET  SIGNAL 
C 

CALL  STDEV(DYO,NCYC,DBAR, SDEV) 

C 

C  DETERMINE  NORMALIZED  KURTOSIS  OF  NET  SIGNAL 
C 

CALL  KURTOS (DYO , NCYC , DBAR , RKURT) 

C 

C  WRITE  RESULTS  TO  FILES 

C 

DO  777,  I -1, NCYC 
WRITE(8 , 700)  TO(I) ,FYO(I) 

777  CONTINUE 
C 

DO  778,  I-l.NCYC 
WRITE(9,700)  TO(I) ,RFYO(I) 

778  CONTINUE 

C 

DO  779,  I-l.NCYC 

WRITE (10, 700)  T0(1) ,DYO(I) 

779  CONTINUE 

C 

DO  780,  I-l.NCYC 

URITE(ll , 700)  TO(I),DSYO(I) 

780  CONTINUE 
C 

WRITE(2 ,200)  TIHE(J), SDEV, RKURT 
200  P0RMAT(F10 . 2 ,  F12 . 6 ,  F12 . 6 ) 

C 

C  CLOSE  LOOP  DEPENDANT  FILES 
C 

CLOSE(  5 ,  IOSTAT-IOERR) 

1F(I0ERR  OT.  0)  GOTO  999 
CLOSE (6 . IQSTAT-IOERR) 

IF(IOERR  .GT.  0)  GOTO  999 
CL0SE(7 /IOSTAT-IOERR) 

IF(IOERR  .GT.  0)  GOTO  999 
OOSE(8 ,  IOSTAT-IOERR) 
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0221  IFUOERR  .GT.  0)  GOTO  999 

0222  CLOSE (9 , IOSTAT-IOERR) 

0223  IF(IOERR  .GT.  0)  GOTO  999 

0224  CLOSE(10, IOSTAT-IOERR) 

0225  IF(IOSRR  .GT.  0)  GOTO  999 

0226  CL0SE(11, IOSTAT-IOERR) 

0227  IF(IOERR  .GT.  0)  GOTO  999 

0228  C 

0229  C  CLOSE  LOOP,  MOVE  TO  NEXT  DATA  POINT 
0230  C 

0231  111  CONTINUE 

0232  C 
0233  C 

0234  PRINT  *,  'FM4  VALUES  STORED  IN  FILE  :\FM4F 

0235  STOP 

0236  999  PRINT  *,  'ERROR  NO.  \I0ERR 

0237  END 

0238  C 
0239  C 

0240  C  END  OF  MAIN  PROGRAM 
0241  C 
0242  C 

0243  C**************^^*****************************^********^**-*******^ 
0244  C 

0245  C  ROUTINE:  IFOUR 

0246  C 

0247  C  THIS  SUBROUTINE  PERFORMS  AN  INVERSE  DISCRETE  FOURIER 

0248  C  TRANSFORM  ON  A  DFT  TRACE  USING  MAGNITUDE  AND  PHASE  DATA. 

0249  C 
0250  C 

0251  c****^^*******^****^^^******^**^**^*************^*^******^ 

0252  C 

0253  SUBROUTINE  IF0UR( F , FMAG , PHAS , T , Y) 

0254  C 

0255  DIMENSION  T(1024) ,Y<1024) 

0256  DIMENSION  F(401)  ,DFTRE(401),DmM(401),FMAG(401).PHAS(401) 

0257  DIMENSION  FQ<20) 

0258  COMPLEX  X(1024) 

0259  C 

0260  C.. DEFINE  DFT  SETTINGS 
0261  C 

0262  NPTS  -  1024 

0263  NPOW  -  10 

0264  NFREQ  -  401 

0265  FSPAN  -  F<401)  *  F(l) 

0266  PERIOD  -  FLOAT (  (NFREQ- 1)  )  /  FSPAN 

0267  BW  -  F(2)  *  F(l) 

0268  DELTA  -  PERIOD  /  FLOAT (  (NPTS-1)  ) 

0269  PRINT  *,  PERIOD 

0270  PRINT  *,  BU 

0271  C 

0272  C  INITIALIZE  TIME,  AND  DFT  ARRAYS 
0273  C 

0274  DO  100  I-l.NPTS 

0275  T(I)  -  DELTA  *  FLOAT (  1*1  ) 
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0276 
0277 
0278 
0279 
0280 
0281 
0282 
0283 
0284 
0285 
0286 
0287 
0288 
0289 
0290 
0291 
0292 
0293 
0294 
0295 
0296 
0297 
0298 
0299 
0300 
0301 
0302 
0303 
0304 
0305 
0306 
0307 
0308 
0309 
0310 
0311 
0312 
0313 
0314 
0315 
0316 
0317 
03 1H 
C3)9 
0320 
0321 
0)22 
0323 
0324 
0325 
0326 
0327 
0328 
0329 
0330 


X(I)  -  (  0. ,0.  ) 

100  CONTINUE 
C 

C  CONVERTING  PHASE  AND  MAGNITUDE  TO  REAL  AND  IMAGINARY  FFT  VALUES 
C 

PI  -  3.1415927 
DO  300,  I-l.NFREQ 
RAD  -  PHAS(I)*PI/180. 

DFTRE(I)  -  FMAG(I)  *  COS(RAD) 

DFTIM(I)  -  FMAG(I)  *  SIN(RAD) 

300  CONTINUE 
C 

C . . PERFORM  INVERSE  DFT  ON  REAL  PART 
C 

I FLAG  -  1 
DO  2200  I-l.NFREQ 

X(I)  -  CMPLX(  DFTRE(I) ,0.  ) 

2200  CONTINUE 

CALL  FFT(  IFLAG.NPOW.X  ) 

DO  2 “00  I-l.NPTS 

Y(I)  -  REAL(  X<I)  ) 

2300  CONTINUE 

C 

C.. PERFORM  INVERSE  DPT  ON  IMAGINARY  PART 
C 

DO  2400  1-1 ,NFREQ 

X(I)  -  CMPLX(  DFTIM(I) ,0.  ) 

2400  CONTINUE 

DO  2500  I-NFREQ+l ,NPTS 
X(I)  -  (  Q.,0.  ) 

2500  CONTINUE 

CALL  FFT(  IFLAG.NPOU.X  ) 

DO  2600  l-l.NPTS 

Y(I)  -  Y(I)  •  A1MAG(  X(I)  ) 

2600  CONTINUE 

C 

RETURN 

END 

C 

C 

C**************+********6******^***w****+***********************a***c 

v> 

C  ROUTINE:  FFT 
C  C 

ca*******#av^*****#*a*isnM»*a**a*Tir*a-*a'*i*’*a****iJr*ar**ifr»*a***iiniira#a*#**ik*»**c 

c 

SUBROUTINE  FFT(  IFLAG,NPOU,X  ) 

DIMENSION  X{1) ,CS(2) ,MSK(13) 

COMPLEX  X, CXCS. HOLD ,XA 
EQUIVALENCE  (CXCS.CS) 

C 

PI  -  3.141  592  7 
NMAX  -  2**NPOW 

2Z  -  2.  *  PI  *  FLOATUFLAG)  /  FLOAT  (NMAX) 

MSK(l)  -  NMAX  /  2 
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0331 

DO 

15  1-2, NPOW 

0332 

MSK(I)  -  MSK(I-l)  /  2 

0333 

15 

CONTINUE 

0334 

NN 

-  NMAX 

0335 

MM 

-  2 

0336 

C 

0337 

C. .LOOP  OVER  NPOW  LAYERS 

0338 

C 

0339 

DO 

45  LAYER-1, NPOW 

0340 

NN  -  NN  /  2 

0341 

NV  -  0 

0342 

DO  40  I-1.MM.2 

0343 

II  -  NN  *  I 

0344 

C 

0345 

C. .exes  ■ 

-  CEXP(  2*PI*NV*FLOAT ( I FLAG ) /FLOAT ( NMAX )  ) 

0346 

C 

0347 

W  -  FLOAT (NW)  *  ZZ 

0348 

CS<1)  -  COS<  W  ) 

0349 

CS(2)  -  SIN(  V  ) 

0350 

C 

0351 

C.. COMPUTE  ELEMENTS  FOR  BOTH  HALVES  OF  EACH  BLOCK 

0352 

C 

0353 

DO  20  J-l.NN 

0354 

II  -  II  4  1 

0355 

IJ  -  II  •  NN 

0356 

XA  -  CTCS  *  X<1I) 

0357 

X(IX)  -  X(IJ)  -  XA 

0358 

X(U>  -  X(IJ)  +  XA 

0359 

20 

CONTINUE 

0360 

C 

0361 

C. .BUMP 

UP  SERIES  BY  2,  COMPUTE  REVERSE  ADDRESS 

0362 

C 

0363 

DO  25  IOC-2, NPOW 

0364 

LL  ~  NV  -  HSK(LOC) 

0365 

IF(  LL  )  30,35,25 

0366 

25 

NW  -  LL 

0367 

30 

NV  -  MSK(LOC)  ♦  NW 

0368 

GOTO  40 

0369 

35 

NW  -  HSK(L0C+1) 

0370 

40 

CONTINUE 

0371 

MM  -  KM  *  2 

0372 

45 

CONTINUE 

0373 

C 

0374 

G\  .DO 

FINAL  REARRANOSSkTF,  ALSO  t&ttXWt  BY  DELTA. 

0375 

C 

0376 

NW  -  0 

0377 

DO  80  I-l.NMAX 

0378 

NW1  -  NW  +  1 

0379 

HOLD  -  X(NWl) 

0380 

1F(  NW1* I  }  60,55,50 

0381 

50 

X<NW1)  -  X{I) 

0382 

55 

X(I)  -  HOLD 

0383 

C 

0384 

C. .BUMP  UP  SERIES  BY  l.  AND  COMPUTE  REVERSE  ADDRESS 

0385 

C 

APPB 

03*6  60  DO  65  LOC-l.HPOW 

0387  LL  -  NW  •  MSK(LOG) 

0388  IF(  LL  )  70,75,65 

0389  65  NW  -  LL 

0390  70  NW  -  MSK(LOC)  +  NW 

0391  GOTO  SO 

0392  75  NW  -  MSK(LOC+l) 

0393  80  CONTINUE 

0394  RETURN 

0395  END 

0396  C 
0397  C 

0398  C********************************************************************C 


0399  C  C 
040Q  C  ROUTINE:  FILTER  C 
0401  C  C 
0402  C  THIS  SUBROUTINE  USES  FILTERS  A,  B,  AND  C  DEFINED  IN  THE  MAIN  C 
0403  C  PROGRAM  TO  FILTER  OUT  ANY  NON -GEAR -RELATED  SIGNALS  FROM  THE  C 
0404  C  ANALYSIS.  (BEARING  NOISE,  SLAVE  GEAR  NOISE,  Etc.)  C 
0405  C  C 
0406  C  C 


0407  C****  ****************************  *************s*******  ******■»  ****v.'***C 

0408  c 

0409  SUBROUTINE  FI LTER ( FREQ , FHAO , FFMAG , A FS , APS , BPS , BPS ,  CFS , CFE ) 

0410  C 

0411  DIMENSION  FREQ(40I) , FHAG(401> , FFHAG(40I) 

0412  C 

0413  DO  100,  1-1.401 

0414  FFKAG(I)  -  FHAG(I) 

0415  IF((FREQ(I)  .CE.  AFS)  .AND.  (FREQ(I)  .LE.  APE))  THEN 

0416  FFMAG(l)  -  0.0 

0417  ENDIF 

0418  1F((FREQ(I)  .GE.  BPS)  .AND.  (FREQ(I)  .LE.  BFE))  THEN 

0419  FFHAG(X)  -  0.0 

0420  ENDIF 

0421  IF((PREQ(I)  CE.  CPS)  .AND.  (FREQ(I>  .US.  CPE))  THEN 

0422  FFMAG(I)  -  0.0 

0423  ENDIF 

0424  100  CONTINUE 

0425  RETURN 

0426  END 

0427  G 

0428  C 

0429  u**  *******************  ******  **********  ******  -  ***i**r*****i»******>;*****c 


0430  C  C 
0431  C  ROUTINE:  FFH4  C 
0432  C  C 
0433  C  THIS  SUBROUTINE  FILTERS  OUT  ALL  AMPLITUDES  EXCEPT :  PRIMARY  C 
0434  C  MESH  FREQUENCY  AND  FIRST  HARMONIC  AND  THEIR  FIRST  CORRESPONDING  C 
0435  C  SIDEBANDS,  AND  PRIMARY  SHAFT  FREQ.  (PER  STEWART'S  1977  PAPER)  C 
0436  C  C 
0437  C  C 


0438  c****************  ***********  ****************************************<• 

0439  C 

0440  SUBROUTINE  FFK4 ( FREQ ,  FFMAG , PRIM , SPD , RMAC) 
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0441 

0442 

0443 

0444 

0445 

0446 

0447 

0448 

0449 

0450 

0451 

0452 

0453 

0454 

0455 

0456 

0457 

0458 

0459 

0'v60 

0461 

0462 

0463 

0464 

0465 

0466 

0467 

0468 

0469 

0470 

0471 

0472 

0473 

0474 

0475 

0476 

0477 

0478 

0479 

0480 

0481 

0482 

0483 

0484 

0485 

0486 

0487 

0488 

0489 

0490 

0491 

0492 

0493 

0494 

0495 


C 

DIMENSION  FREQ(401) . FFMAG(401) ,SMAG(401) 

C 

HARM  -  PRIM  *  2. 

LPRIM  -  PRIM  -  1.5  *  SPD 
HPRIM  »  PRIM  +  1.5  *  SPD 
LHARM  -  HARM  -  1.5  *  SPD 
HHARM  -  HARM  +  1.5  *  SPD 
I, SHAFT  -  0. 5  *  SPD 
HSHAFT ■-  1.5  *  SPD 
C 

DO  100,  1-1,401 
RMAG(I)  -  0.0 

IF( (FREQ(I)  .GT.  LPRIM)  .AND.  (FREQ(I)  .LT.  HPRIM))  THEN 
RMAG(I)  -  FFMAG(I) 

END  IF 

IF( (FREQ(I)  .GT.  LHARM)  .AND.  (FREQ(I)  .LT.  HHARM))  THEN 
RMAG(I)  -  FFMAG(I) 

ENDIF 

IF( (FREQ(I)  .GT.  LSHAFT)  .AND.  (FREQ(I)  .LT.  HSHAFT))  THEN 
RMAG(I)  -  FFMAG(I) 

ENDIF 

100  CONTINUE 
RETURN 
END 
C 
C 

C********************************************************************C 


C  c 

C  ROUTINE:  STDEV  C 
C  C 
C  THIS  SUBROUTINE  COMPUTES  THE  STANDARD  DEVIATION  OF  THE  C 
C  DIFFERENCE  LIME  SIGNAL  (ORIGINAL  -  REGULAR  COMPONENTS).  C 
C  C 
C  C 


C********************************************************************(j 

C 

SUBROUTINE  STDEV <DYO,N,DBAR,SD) 

C 

DIMENSION  DY0(220) 

ASUM  -  0.0 
DO  100,  I»1,N 

ASUM  -  (DYO(I)  -  DBAR)**2 .  +  ASUM 
100  CONTINUE 

SO  -  ( ASUM/FLOAT ( N) ) ** , 5 
RETURN 
END 
C 
C 

C*********** ************* **************************************  ******C 


C  C 
C  ROUTINE:  KURTOS  C 
C  C 
C  THIS  SUBROUTINE  COMPUTES  THE  NORMALIZED  KURTOSIS  OF  THE  C 
C  DIFFERENCE  TIME  SIGNAL  v ORIGINAL  -  REULAR  COMPONENTS)  C 
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0496  C 
0497  C 

0498  C**w***********************************-Jr*****************************c 


0499 

C 

0500 

SUBROUTINE  KURTOS (DYQ , N , DBAR , RKT ) 

0501 

G 

0502 

DIMENSION  DYO(220) 

0503 

C 

0504 

TSUM  -  0.0 

0505 

DO  100,  I-l.N 

0506 

TSUM  -  TSUM  +  (DYO(I)  -  DBAR)**4. 

0507 

100 

CONTINUE 

0508 

C 

0509 

BSUM  -  0.0 

0510 

DO  200,  I-1,N 

0511 

BSUM  -  BSUM  +  (DYO(I)  -  DBAR)**2. 

0512 

200 

CONTINUE 

0513 

C 

0514 

RKT  -  { FLOAT (N)  *  TSUM) /BSUM**2 . 

0515 

RETURN 

0516 

END 
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APPENDIX  C  -  SOURCE  CODE  FOR  HUB. FOR 


OOOl 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 

0011 

0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 

0021 

0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 

0051 

0052 

0053 

0054 

0055 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 


c 


c 

PROGRAM:  HILR . FOR  REVISION:  A  C 

REVISION  DATE:  05-12-89  C 
PROGRAMMER:  J.J.  Zakrajsek  DATE  CREATED:  05-03-89  C 

C 

THIS  PROGRAM  USES  THE  HILBERT  TRANSFORM  THEORY  TO  CONVERT  C 
THE  TIME  SIGNAL  TO  A  COMPLEX  FUNCTION  TO  DEMODULATE  THE  C 
SIDEBANDS  FROM  THE  CARRIER  FREQUENCY.  THIS  METHOD  IS  BASED  C 
ON  THE  1986  PAPER  BY  P.  D.  McFADDEN ,  AND  THE  1988  PAPER  BY  C 
D.  A.  WALLACE  AND  M.  K.  DARLOW.  C 

C 

REV  A:  ADD  VARIABLE  RAMP  FREQUENCY  ANALYTIC  PHASE  VALUE  C 
CORRECTION  FUNCTION.  C 

C 


C 

DIMENSION  TIME(50) ,AMAGN(220) ,APHAS(220) ,PPHAS(220) ,CPHAS(220) 
DIMENSION  FREQ(401) ,FMAGRMS(401) ,FMAG(401) ,PHAS(401) ,FFMAG(401) 
DIMENSION  T0(1024) ,YACT(1024) ,YHIL(1024) ,POINT(60) 

CHARACTER*22  LIST1 , LI ST 3 , AMPLF , PHASEF 
CHARACTER*22  AAMPLF , APHASF 
COMPLEX  X(1024) 

C 

C  ENTER  GENERAL  PARAMETERS  AND  FILE  NAMES 
C 

PRINT  *,  'ENTER  NUMBER  OF  DATA  POINTS  (INTEGER)' 

READ(*,*)  K 

PRINT  *,  'ENTER  CARRIER  FREQUENCY  (PRIMARY  OR  FIRST  HARMONIC)' 
READ(*,*)  PRIM 

PRINT  *,  'ENTER  ROTATIONAL  SPEED  IN  REV/SEC  (Hz)' 

READ(*,*)  SPD 

PRINT  *,  'ENTER  NUMBER  OF  TIME  DATA  POINTS  PER  CYCLE' 

PRINT  *,  '(  INTEGER,  220  MAXIMUM  )' 

READ  (*,*)  NCYC 

C 

C  DEFINITION  OF  FREQUENCY  BAND  TO  BE  USED  IN  ANALYSIS 
C 

PRINT  *,  'THIS  ANALYSIS  IS  BASED  ON  A  BAND  OF  FREQUENCY  DATA' 
PRINT  *,  'SURROUNDING  THE  MESH  FREQUENCY  CHOSEN  ABOVE.  THE  ' 
PRINT  *,  'NUMBER  OF  SIDEBANDS  ON  EACH  SIDE  OF  THE  FREQUENCY' 
PRINT  *,  'OF  INTEREST  ARE  CHOSEN  BY  THE  USER.  DATA  OUTSIDE  THIS' 
PRINT  *,  'BAND,  ALONG  WITH  THE  FREQUENCY  CHOSEN  ABOVE,  ARE' 

PRINT  *,  'FILTERED  OUT  OF  THE  ANALYSIS.' 

PRINT  *,  ' 

C 

PRINT  *,  'ENTER  NUMBER  OF  SIDEBANDS  TO  BE  PASSED  BELOW' 

PRINT  *,  'CHOSEN  FREQUENCY  (  REAL  NUMBER  XX.)' 

READ(*,*)  FILT3 

PRINT  *,  'ENTER  NUMBER  OF  SIDEBANDS  TO  BE  PASSED  ABOVE' 

PRINT  *,  'CHOSEN  FREQUENCY  (REAL  NUMBER  XX.)' 

READ(*,*)  FILTE 
C 

C  INPUT  PRIMARY  FILE  NAMES  FILES 
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✓ 


0056  C 

0057  PRINT  *,  'ENTER  TXXLISTl.FIL  INPUT  FILES  FILENAME' 

0058  READ(*,'(A22)')  LIST1 

0059  PRINT  *,  'ENTER  TXXLIST3.FIL  OUTPUT  FILES  FILENAME' 

0060  READ(*, ' (A22) ' )  LIST3 

0061  C 

0062  C  OPEN  PRIMARY  FILES 
0063  C 

0064  OPEN ( 3 , FILE-LIST1 , IOSTAT-IOERR) 

0065  IF(IOERR  .GT.  0)  GOTO  999 

0066  0PEN(4,FILE-LIST3, IOSTAT-IOERR) 

0067  IF(IOERR  .GT.  0)  GOTO  999 

0068  C 

0069  C  STARTING  POINT  OF  ANALYSIS  LOOP  FOR  EACH  DATA  POINT 
0070  C 

0071  DO  111,  J-l.K 

0072  PRINT*,  'PROCESSING  DATA  POINT  NO.  :',J 

0073  READ(3 , 300)  AMPLF.PHASEF 

0074  300  FORMAT(A22,A22) 

0075  READ(4,400)  AAMPLF, APHASF 

0076  400  FORMAT ( A2 2 , A2 2 ) 

0077  C 

0078  C  OPEN  LOOP  DEPENDANT  FILES 
0079  C 

0080  OPEN ( 5, FILE-AMPLF, IOSTAT-IOERR) 

0081  IF(IOERR  .GT.  0)  GOTO  999 

0082  OPEN (6 , FILE-PHASEF, IOSTAT-IOERR) 

0083  IF(IOERR  .GT.  0)  GOTO  999 

0084  OPEN ( 8, FILE-AAMPLF, IOSTAT-IOERR) 

0085  IF(IOERR  .GT.  0)  GOTO  999 

0086  0PEN(9 , FILE-APHASF, IOSTAT-IOERR) 

0087  IF(IOERR  .GT.  0)  GOTO  999 

0088  C 

0089  C  READ  IN  PHASE  AND  FREQUENCY  FILES 
0090  C 

0091  DO  222,  1-1,401 

0092  READ(5 , 500)  FREQ(I) ,FMAGRMS(I) 

0093  500  F0RMAT(4X , F12 . 3 , 20X , F16 . 8 ) 

0094  FMAG(I)  -  FMAGRMS(I)  *  (2. **.5) 

0095  READ(6,600)  PHAS(I) 

0096  600  FORMAT (20X.F12. 3) 

0097  222  CONTINUE 

0098  C 

0099  C  PERFORM  BANDPASS  FILTERING 
0100  C 

0101  CALL  FILTER(FREQ , FMAG , FFMAG , FILTS , FILTE , PRIM , SPD) 

0102  C 

0103  C  PERFORM  IFFT  ON  FILTERED  FREQUENCY  COMPONENTS  OF  ACTUAL  SIGNAL 
0104  C 

0105  IFL  -  1 

0106  CALL  IFOUR(FREQ, FFMAG, PHAS, TO, YACT, IFL) 

010^  C 

0108  C  OBTAIN  HILBERT  TRANSFORM  OF  ACTUAL  SIGNAL 
0109  C 

0110  C  PERFORM  IFFT  ON  FREQUENCY  COMPONENTS  OF  HILBERT  TRANSFORM  OF  SIGNAL 
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0111 

0112 

0113 

0114 

0115 

0116 

0117 

0118 

0119 

0120 

0121 

0122 

0123 

0124 

0125 

0126 

0127 

0128 

0129 

0130 

0131 

0132 

0133 

0134 

0135 

0136 

0137 

0138 

0139 

0140 

0141 

0142 

0143 

0144 

0145 

0146 

0147 

0148 

0149 

0150 

0151 

0152 

0153 

0154 

0155 

0156 

0157 

0158 

0159 

0160 

0161 

0162 

0163 

0164 

0165 


C 

IFL  -  2 

CALL  IFOUR ( FREQ , FFMAG , PHAS , TO , YHIL , IFL) 

C 

C  DETERMINATION  OF  ANALYTIC  FUNCTION  MAGNITUDE  AND  PHASE 
C 

PI  -  3.1415927 
DO  444,  I-l.NCYC 

AMAGN(I)  -  (YACT(I)**2 .  +  YHIL(I)**2 .)**(. 5) 

APHAS(I)  -  (ATAN2(YHIL(I) ,YACT(I) ) )*(180 ./PI) 

444  CONTINUE 
C 

C  CORRECT  PHASE  VALUES  FOR  NOMINAL  GEAR  MESH  FREQUENCY 
C 

C  ADD  180  DEGREES  TO  ALL  PHASE  VALUES 
C 

DO  210,  I-l.NCYC 
PPHAS(I)  -  APHAS(I)  +  180. 

210  CONTINUE 
C 

C  DETERMINE  FREQUENCY  OF  RAMP  FUNCTION 
C 

NCT  -  0 

LIM  -  NCYC  -  1 
DO  211,  I-l.LIM 
N  -  I  +  1 

XCHK  -  PPHAS(I)  -  p£hAS(N) 

IF(XCHK  .GT.  200.)  THEN 
NCT  -  NCT  +  1 

POINT (NCT)  -  (TO(I)  +  TO(N))/2. 

END  IF 

211  CONTINUE 

SFREQ  -  FLOAT (NCT-l)/( POINT (NCT)  -  POINT(l)) 

PRINT  *,  SFREQ 
C 

C  DETERMINE  INITIAL  OFFSET  OF  RAMP  FUNCTION 
C 

I  -  0 

212  CONTINUE 
I  -  I  +  1 
N  -  I  +  1 

XCHK  -  PPHAS(I)  -  PPHAS(N) 

IF (XCHK  .GT.  200.)  GOTO  213 
GOTO  212 

213  PERIOD  -  1 . 0/SFREQ 

OFFSET  -  PERIOD  -  (TO(I)  +  TO(N))/2.0 
PRINT  *,  OFFSET 
C 

C  CORRECT  PHASE  ANGLE  VARIATION  FOR  NOMINAL  TOOTH  MESH  FREQUENCY 
C 

DO  214,  I-l.NCYC 
IF(I  .EQ.  1)  THEN 
FC  -  1.0 
ENDIF 

IF(I  .EQ.  1)  GOTO  379 
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0166  N  -  I  -  1 

0167  XCHK  -  PPHAS(N)  -  PPHAS(I) 

0168  IF(XCHK  .GT.  200.)  THEN 

0169  M  -  I  -  1 

0170  377  M  -  M  +  1 

0171  MN  -  M  +  1 

0172  CCHK  -  PPHAS(M)  -  PPHAS(MN) 

0173  IF (CCHK  .GT.  200.)  GOTO  378 

0174  GOTO  377 

0175  378  LL  -  I  -  1 

0176  ULIM  -  (TO(MN)  +  T0(M))/2.0 

0177  BLIM  -  (TO(I)  +  TO(LL))/2.0 

0178  C  PRINT  *,  ULIM 

0179  C  PRINT  *,  BLIM 

0180  CFREQ  -  1.0/ (ULIM  -  BLIM) 

0181  FC  -  CFREQ/SFREQ 

0182  OFFSET  -  -1.0*BLIM 

0183  ENDIF 

0184  379  CONTINUE 

0185  CPHAS(I)  -  PPHAS(I)  -  (360 . )*SFREQ*FC*(TO(I)  +  OFFSET) 

0186  214  CONTINUE 

0187  C 

0188  C  WRITE  RESULTS  TO  FILES 
0189  C 

0190  DO  777,  I-l.NCYC 

0191  WRITE(8 , 700)  TO(I) , AMAGN(I) 

0192  700  F0RMAT(F16.8,4X,F16.8) 

0193  777  CONTINUE 

0194  C 

0195  DO  778,  I-l.NCYC 

0196  WRITE(9,700)  TO(I) ,CPHAS(I) 

0197  778  CONTINUE 

0198  C 

0199  C  CLOSE  LOOP  DEPENDANT  FILES 
0200  C 

0201  CLOSE(5 , IOSTAT-IOERR) 

0202  IF(IOERR  .GT.  0)  GOTO  999 

0203  CL0SE(6 , IOSTAT-IOERR) 

0204  IF(IOERR  .GT.  0)  GOTO  999 

0205  CLOSE (8, IOSTAT-IOERR) 

0206  IF(I0ERR  .GT.  0)  GOTO  999 

0207  CLOSE(9, IOSTAT-IOERR) 

0208  IF(I0ERR  .GT.  0)  GOTO  999 

0209  C 

0210  C  CLOSE  LOOP,  MOVE  TO  NEXT  DATA  POINT 
0211  C 

0212  111  CONTINUE 

0213  C 
0214  C 

0215  PRINT  *,  'VALUES  STORED  IN  FILES  ON  DISK  :A' 

0216  STOP 

0217  999  PRINT  *,  'ERROR  NO.  ' , I0ERR 

0218  END 

0219  C 
0220  C 
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0221 

0222 

0223 

0224 

0225 

0226 

0227 

0228 

0229 

0230 

0231 

0232 

0233 

0234 

0235 

0236 

0237 

0238 

0239 

0240 

0241 

0242 

0243 

0244 

0245 

0246 

0247 

0248 

0249 

0250 

0251 

0252 

0253 

0254 

0255 

0256 

0257 

0258 

0259 

0260 

0261 

0262 

0263 

0264 

0265 

0266 

0267 

0268 

0269 

0270 

0271 

0272 

0273 

0274 

0275 


C  END  OF  MAIN  PROGRAM 

C 

C 


C  C 
C  ROUTINE:  I FOUR  C 
C  C 
C  THIS  SUBROUTINE  PERFORMS  AN  INVERSE  DISCRETE  FOURIER  C 
C  TRANSFORM  ON  A  DFT  TRACE  USING  MAGNITUDE  AND  PHASE  DATA.  C 
C  C 
C  G 


C**^^*^^***'iHf^*SHf**VM:****^ihHNriNr^*'********************'********C 


C 

SUBROUTINE  I FOUR ( F , FMAG , PHAS , T , Y , I FL) 

C 

DIMENSION  T(1024) ,Y(1024) 

DIMENSION  F(401) ,DFTRE(401) ,DFTIM(401) ,FMAG(401) ,PHAS(401) 
DIMENSION  FQ{20) 

COMPLEX  X(1024) 

C 


C. .DEFINE  DFT  SETTINGS 
C 

NPTS  -  1024 

NPOW  -  10 

NFREQ  -  401 

FSPAN  -  F(401)  -  F(l) 

PERIOD  -  FLOAT (  (NFREQ- 1)  )  /  FSPAN 
BW  -  F(2)  -  F(l) 

DELTA  -  PERIOD  /  FLOAT (  (NPTS-1)  ) 

PRINT  *,  PERIOD 
PRINT  *,  BW 

C 

C  INITIALIZE  TIME,  AND  DFT  ARRAYS 
C 

DO  100  I-l.NPTS 

T(I)  -  DELTA  *  FLOAT (  1-1  ) 

X(I)  -  (  0..0.  ) 

100  CONTINUE 
C 

C  CONVERTING  PHASE  AND  MAGNITUDE  TO  REAL  AND  IMAGINARY  FFT  VALUES 
C 

PI  »  3.1415927 
C 

G  ACTUAL  SIGNAL  ALGORYTHM  (IFL  -  1) 

C 

IF(IFL  .EQ.  1)  THEN 
DO  300,  1-1, NFREQ 
RAD  -  PHAS(I)*PI/180. 

DFTRE(I)  -  FMAG(I)  *  COS(RAD) 

DFTIM(I)  -  FMAG(I)  *  SIN(RAD) 

300  CONTINUE 
ENDIF 
C 

C  HILBERT  TRANSFORM  ALGORYTHM  (IFL  -  2) 

C 
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0276 

0277 

0278 

0279 

0280 

0281 

0282 

0283 

0284 

0285 

0286 

0287 

0288 

0289 

0290 

0291 

0292 

0293 

0294 

0295 

0296 

0297 

0298 

0299 

0300 

0301 

0302 

0303 

0304 

0305 

0306 

0307 

0308 

0309 

0310 

0311 

0312 

0313 

0314 

0315 

0316 

0317 

0318 

0319 

0320 

0321 

0322 

0323 

0324 

0325 

0326 

0327 

0328 

0329 

0330 


IF(IFL  .EQ.  2)  THEN 
DO  303,  I-l.NFREQ 
RAD  -  PHAS(I)*PI/180. 

DFTRE(I)  -  FMAG(I)*SIN(RAD) 
DFTIM(I)  -  (-1.0)*FMAG(I)*C0S(RAD) 
303  CONTINUE 
ENDIF 


C. .PERFORM  INVERSE  DFT  ON  REAL  PART 
C 

I FLAG  -  1 
DO  2200 -I-l.NFREQ 

X(I)  -  CMP LX (  DFTRE(I) ,0.  ) 
2200  CONTINUE 

CALL  FFT(  IFLAG.NPOW.X  ) 

DO  2300  I-l.NPTS 

Y(I)  -  REAL(  X(I)  ) 

2300  CONTINUE 


C. .PERFORM  INVERSE  DFT  ON  IMAGINARY  PART 
C 

DO  2400  I-l.NFREQ 

X(I)  -  CMPLX(  DFTIM(I) ,0.  ) 
2400  CONTINUE 

DO  2500  I-NFREQ+1 , NPTS 
X(I)  -  (  0.,0.  ) 

2500  CONTINUE 

CALL  FFT<  IFLAG,NPOW.X  ) 

DO  2600  I-l.NPTS 

Y(I)  -  Y(I)  -  AIMAG(  X(I)  ) 
2600  CONTINUE 

C 

RETURN 

END 

C 

C 


C 

C 

C 


ROUTINE:  FFT 


SUBROUTINE  FFT(  IFLAG.NPOW.X  ) 

DIMENSION  X(l) ,CS(2) ,MSK(13) 

COMPLEX  X, CXCS, HOLD, XA 
EQUIVALENCE  (CXCS.CS) 

C 

PI  -  3,141  592  7 
NMAX  -  2**NPOW 

ZZ  -  2.  *  PI  *  FLOAT(IFLAG)  /  FLOAT(NMAX) 
MSK(l)  -  NMAX  /  2 
DO  15  I-2.NPOW 

MSK(I)  -  MSK(I-l)  /  2 
15  CONTINUE 

NN  -  NMAX 
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0331 

MM 

-  2 

0332 

C 

0333 

C. .LOOP  OVER  NPOW  LAYERS 

0334 

C 

0335 

DO 

45  LAYER-1, NPOW 

0336 

NN  -  NN  /  2 

0337 

NW  -  0 

0338 

DO  40  I-1.MM.2 

0339 

II  -  NN  *  I 

0340 

C 

0341 

C. .GXCS  ■ 

-  CEXP(  2*PI*NW*FL0AT(IFLAG)/FL0AT(NMAX)  ) 

0342 

G 

0343 

W  -  FLOAT (NW)  *  ZZ 

0344 

CS(1)  -  C0S(  W  ) 

0345 

CS(2)  -  SIN(  W  ) 

0346 

C 

0347 

C. .COMPUTE  ELEMENTS  FOR  BOTH  HALVES  OF  EACH  BLOCK 

0348 

C 

0349 

DO  20  J-l.NN 

0350 

II  -  II  +  1 

0351 

IJ  -  II  -  NN 

0352 

XA  -  CXCS  *  X(II) 

0353 

X(II)  -  X(IJ)  •  XA 

0354 

X(IJ)  -  X(IJ)  +  XA 

0355 

20 

CONTINUE 

0356 

C 

0357 

C. .BUMP  1 

UP  SERIES  BY  2,  COMPUTE  REVERSE  ADDRESS 

0358 

C 

0359 

DO  25  LOC-2 , NPOW 

0360 

LL  -  NW  -  MSK(LOC) 

0361 

IF(  LL  )  30,35,25 

0362 

25 

NW  -  LL 

0363 

30 

NW  -  MSK(LOC)  4-  NW 

0364 

GOTO  40 

0365 

35 

NW  -  KSK(LOC+l) 

0366 

40 

CONTINUE 

0367 

MM  -  MM  *  2 

0368 

45 

CONTINUE 

0369 

C 

0370 

C.  .DO 

FINAL  REARRANGEMENT,  ALSO  MULTIPLY  BY  DELTA. 

0371 

C 

0372 

NW 

-  0 

0373 

DO 

80  I-l.NMAX 

0374 

NW1  -  NW  +  1 

0375 

HOLD  -  X(NW1) 

0376 

IF(  NW1-I  )  60,55,50 

0377 

50 

X(NWl)  -  X(I) 

0378 

55 

X(I)  “  HOLD 

0379 

C 

0380 

C. .BUMP  UP  SERIES  BY  l,  AND  COMPUTE  REVERSE  ADDRESS 

0381 

C 

0382 

60 

DO  65  LOC-1 , NPOW 

0383 

LL  -  NW  -  MSK(LOC) 

0384 

IF(  LL  )  70,75,65 

0385 

65 

NW  -  LL 
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0386 

0387 

0388 

0389 

0390 

0391 

0392 

0393 

0394 

0395 

0396 

0397 

0398 

0399 

0400 

0401 

0402. 

0403 

0404 

0405 

0406 

0407 

0408 

0409 

0410 

0411 

0412 

0413 

0414 

0415 

0416 

0417 

0418 

0419 

0420 

0421 

0422 

0423 

0424 


70  NW  -  MSK(LOC)  +  NW 

GOTO  80 

75  NW  -  MSK(L0C-+1) 

80  CONTINUE 

RETURN 
END 
C 
C 


C 

C 

C 

C 

C 

C 

C 

C 


ROUTINE:  FILTER 

THIS  SUBROUTINE  FILTERS  OUT  ALL  FREQUENCY  COMPONENTS 
OUTSIDE  THE  FREQUENCY  BAND  CHOSEN  IN  THE  MAIN  PROGRAM, 
THE  CARRIER  FREQUENCY  IS  ALSO  FILTERED  OUT. 


C 

C 

c 

c 

c 

c 

c 

c 


c 

SUBROUTINE  FILTER ( FREQ , FMAG , FFMAG , FS , FE , PRIM , SPD) 

C 

DIMENSION  FREQ<401) . FMAG(4Q1) , FFMAG (401) 

C 

SBL  -  PRIM  -  (FS  +  .5)*SPD 
PL  -  PRIM  -  ( . 5*SPD) 

PU  -  PRIM  +  ( . 5*SPD) 

SBU  -  PRIM  4-  (FE  +  . 5)*SPD 
C 

DO  100,  1-1,401 
FFMAG(I)  -  0.0 

IF( (FREQ(I)  .GE.  SBL)  .AND.  (FREQ(I)  .LE.  PL))  THEN 
FFMAG ( I )  -  FMAG(I) 

END  IF 

IF( (FREQ(I)  .GE.  PU)  .AND.  (FREQ(I)  -LE.  SBU))  THEN 
FFMAG(I)  -  FMAG(I) 

ENDIF 

100  CONTINUE 
RETURN 
END 
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APPENDIX  D  -  SOURCE  CODE  FOR  PARAM.FOR 

OOOl  c********************************************************************c 


0002  C  c 

0003  C  C 

0004  C  PROGRAM:  PARAM.FOR  CREATED:  05-05-89  REVISED:  XX-XX-XXXX  C 
0005  C  C 

0006  C  PROGRAMMER:  J.J.  ZAKRAJSEK  C 

0007  C  C 

0008  C  THIS  PROGRAM  COMPUTES  SEVERAL  FAULT  DETECTION  PARAMETERS  C 

0009  C  BASED  ON  THE  A-B-C  FILTERED  SIGNAL  CALCULATED  IN  THE  PROGRAM  C 

0010  C  TFM4 . FOR.  THE  BASIS  OF  THE  PARAMETERS  ARE  SUMMARIZED  BELOW.  C 

0011  C  C 

0012  C  CREST  FACTOR  (CF)  -  PEAK  LEVEL  /  TOTAL  RMS  LEVEL  C 

0013  C  C 

0014  C  SIDEBAND  LEVEL  FACTOR  (SBLF)  -  FIRST  ORDER  SIDEBANDS  (PRIMARY)  C 

0015  C  TOTAL  RMS  LEVEL 

0016  C 

0017  C  WEAR  ENERGY  RATIO  -  RMS  LEVEL  OF  NON-PERIODIC  PART  OF  SIGNAL 

0018  C  (WER)  /  RMS  LEVEL  OF  PERIODIC  PART  OF  SIGNAL 

0019  C 

0020  C 


0021  C********************************************************************C 

0022  C 

0023  DIMENSION  TIME(50) , FSIG(220) ,RSIG(220) ,DSIG(220) 

'024  DIMENSION  FREQ(500),  AMPL(500) 

0025  CHARACTER*22  FMOF, PARAMF, LIST2 , LIST1 , AHPLF 

0026  CHARACTER* 18  FT1MEF.RTIMEF.DTIMEF 

002  7  C 

0028  C  ENTER  INITIAL  DATA  AND  FILE  NAMES 
0029  C 

0010  PRINT  *,  'ENTER  NUMBER  OF  DATA  POINTS* 

0031  READ<*.*)  K 

0032  PRINT  *,  'ENTER  MESH  FREQ  (XJXXX.)* 

0033  READ(*,*)  PRIM 

0034  PRINT  *.  'ENTER  ROTATIONAL  SPEED  (REV/SEC)' 

0035  READ(*,*)  SPD 

0036  C 

0037  PRINT  *.  'ENTER  OUTPUT  DATA  FILE  NAME  (\TXXPARAH.DAT)* 

0033  READ  (*,*<A22)*)  PARAMF 

0039  OPEN ( 7 , FILE- PARAMF, IOSTAT-IOERR) 

0040  IF  (IOERR  .GT.  0)  GOTO  999 

0041  PRINT  *,  'ENTER  FILE  NAMES  ARRAY  FILE  (\TXXLISTl . FIL) * 

0042  READ(*, ' (A22) ')  LIST1 

0043  OPEN (6 . FILE-LIST1 . IOSTAT-IOERR) 

0044  I F( IOERR  .GT.  0)  GOTO  999 

0045  PRINT  *,  'ENTER  FILE  NAMES  ARRAY  FILE  (\TXXLIST2 . FIL) * 

0046  READ(*. * (A22) * )  LIST2 

004  7  OPF.N<8 ,  FI LE-LIST2  ,  IOSTAT-IOERR) 

0048  IF( IOERR  .GT.  0)  GOTO  999 

0049  PRINT  *,  'ENTER  PREVIOUS  TFMO  FILE  (\TXXFM0.DAT) * 

0050  READ(*,’(A22)')  FMOF 

0051  0PEN(9,FILE-FM0F, IOSTAT-IOERR) 

0052  I F( IOERR  .GT.  0)  GOTO  999 

0053  PRINT  *,  'ENTER  NUMBER  OF  TIME  DATA  POINTS  PER  CYCLE  (INTEGER)' 

0054  READ(*,*)  NCYC 

0055  C 
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0056  C  START  TIME  LOOPING,  PROCESSING  DATA  AT  EACH  TIME  INTERVAL 
0057  C 

0058  DO  300,  J-1,K 

0059  PRINT  *,  'ACCESSING  AND  PROCESSING  DATA  FILE  NO: \J 

0060  READ(6,600)  AMPLF 

0061  600  FORMAT (A20) 

0062  READ(8 , 800)  FTIMEF , RTIMEF , DTIMEF 

0063  800  F0RMAT(A18 , A18 , A18) 

0064  READ(9 , 900)  TIME(J) 

0065  900  FORMAT (F8. 2) 

0066  0PEN(4, FILE-FTIMEF, IOSTAT-IOERR) 

0067  IF  (IOERR  .GT.  0)  GOTO  999 

0068  CPEN(3 .FILE-RTIMEF, IOSTAT-IOERR) 

0069  IF  (IOERR  .GT.  0)  GOTO  999 

0070  OPEN(2 , FILE-DTIMEF, IOSTAT-IOERR) 

0C71  IF  (IOERR  .GT.  0)  GOTO  999 

0072  OPEN(5,FILE-AMPLF, IOSTAT-IOERR) 

0073  IF  (IOERR  .GT.  0)  GOTO  999  . 

0074  C 

0075  C  READ  IN  AND  FIND  THE  MAX  VALUE  OF  THE  FILTERED  TIME  RECORD 
0076  C 

0077  PKL  -  0.0 

0078  DO  100,  I-l.NCYC 

0079  READ(4 ,400)  FSIG(I) 

0080  400  F0RMAT(20X.F16.8) 

0081  AVALU  -  ABS(FS1G(I>) 

0082  IF<  AVALU  .GT.  PKL  )  THEN 

0083  PKL  -  AVALU 

0084  ENDIF 

0085  100  CONTINUE 

0086  C 

0087  C  FIND  FIRST  ORDER  SIDEBANDS  LEVEL 
0088  C 

0089  DO  301,  1-1,401 

0090  READ (5, 500)  FREQ(I) ,AKPL(I) 

0091  500  FORHAT(4X.F12.3.20X,F16.8) 

0092  301  CONTINUE 

0093  C 

0094  CALL  FOSB(FREQ,AMPL,PRIN,SPD.$BL) 

0093  C 

0096  C  FIND  TOTAL  RNS  LEVEL  OF  SIGNAL 
0097  C 

0098  CALL  STDEV ( FS IG , NCYC , FSD ) 

0099  C 

0100  C  FIND  RMS  LEVEL  OF  PERIODIC  COMPONENTS 
0101  C 

0102  DO  120,  I -1, NCYC 

0103  READ(3 ,400)  RSIG(I) 

0104  120  CONTINUE 

0105  C 

0106  CALL  STDEV(RSIG,NCYC,RSD) 

0107  C 

0108  C  FIND  RMS  LEVEL  OF  DIFFERENCE  SIGNAL  (SIGNAL  -  PERIODIC  COMPONENTS) 
0109  C 

0110  DO  140,  I-l.NCYC 
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0111 

0112 

0113 

0114 

0115 

0116 

0117 

0118 

011S 

0120 

0121 

0122 

0123 

0124 

0125 

0126 

0127 

0128 

0129 

0130 

0131 

0132 

0133 

0134 

0135 

0136 

0137 

0138 

0139 

0140 

0141 

0142 

0143 

0144 

0145 

0146 

0147 

0148 

0149 

0150 

0151 

0152 

0153 

0154 

0155 

0156 

0157 

0158 

0159 

0160 

0161 

0162 

0163 

0164 

0165 


READ(2,400)  DSIG(I) 

140  CONTINUE 
C 

CALL  STOEV(DSIG  ,NCYC. ,  DSD) 

C 

C  COMPUTE  CREST  FACTOR 
C 

CF  -  PKL/FSD 
C 

C  FIND  SIDEBAND  LEVEL  FACTOR 
C 

SBLF  -  SBL/FSD 
C 

C  FIND  WEAR  ENERGY  RATIO 
C 

WER  »  DSD/RSD 
C 

C  CLOSE  LOOP  DEPENDANT  FILE,  AND  WRITE  CALCULTED  DATA  TO  FILE 
C 

CL0SE(4 , IOSTAT-IOERR) 

IF  (IOERR  .GY.  0)  GOTO  999 
CLOSE ( 2 , IOSTAT-IOERR ) 

IF  (IOERR  .GT.  0)  GOTO  999 
CLOSE ( 3 , IOSTAT-IOERR) 

IF  (IOERR  .GT.  0)  GOTO  999 
CLOSE( 5. IOSTAT-IOERR) 

IF  (IOERR  .GT.  0)  GOTO  999 
WRITE  (7,700)  TIHE(J),CF, SBLF, WER 
700  FORMAT ( F8 . 2 , F10 . 4 , F10 . 4 , FIG .4) 

300  CONTINUE 
C 

PRINT  *,  ‘PARAMETER  VALUES  STORED  IN  FILE; ' , PARAMF 
STOP 

999  PRINT*.  ‘ERROR  NO, ‘.IOERR 
END 
C 
C 

C  ******  END  OF  MAIN  PROGRAM  ******* 

C 

C 


C 

SUBROUTINE  FOSB(FREQ , AHPL.PRIM . SPD . SBL) 
G 

DIMENSION  FREQ(500) ,AHPL(500) 

BSBS  -  PRIM  *  1.5*SPD 
BSBE  -  PRIM  *  . 5*$PD 
USBS  -  PRIH  +  . 5*$PD 
USBE  -  PRIM  ♦  l.$*SPD 
BSB  -  0.0 
USB  -  0.0 
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0166 
0167 
0168 
0169 
0170 
0171 
0172 
0173 
0174 
0175 
0176 
0177 
0178 
0179 
0180 
0181 
0182 
0183 
0184 
0185 
0186 
0187 
0188 
0189 
0190 
0191 
0192 
0193 
0196 
0195 
0196 
0197 
01 9S 
0199 
0200 


DO  200,  1-1,401 

IF( (FREQ(I) .GT. BSBS)  .AND.  (FREQ(I) .LT.BSBE))  THEN 
IF(AMPL(I)  .GT.  BSB)  BSB  -  AMPL(I) 

END  IF 

IF((FREQ(I)  .GT.USBS)  .AND.  ( FREQ ( I )  . LT . IIS8E-) )  THEN 
IF(AMPI*(I)  .GT.  USB)  USB  -  AMPL(I) 

ENDIF 

200  CONTINUE 

SBL  -  (BSB  +  USB)/2.0 
RETURN 
END 
C 
C 

C 

c  mmim:  stdsv  (standard  deviation  of  time  signal) 

c  .. 

c 

SUBROUTINE  STDEV(DYO,N,SD> 

C 


c 


DIMENSION  0*0(220) 

BSUM  -  0.0 
DO  100,  1-1, N 
DSUH  -  DSUH  +  DVO(I) 

100  CONTINUE 

DBAR  -  DSUH/FLOAT(N) 

ASUM  -  0.0 
DO  200,  I-l.N 

ASUM  -  (DYO(I)  -  DEAR)** 2 .  +  ASUM 
200  CONTINUE 

SD  -  <ASUH/FL0AT{N))**.5 

RETURN 

END 
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APPENDIX  E  -  NOMENCLATURE 


A 

A(t) 

AN(f) 

D(t) 

d 

d 

FCm) 
f  <n> 
F(m) 

f(n) 

F(w) 

f 

f(min) 

FOSL 

Gxy(f> 

Gxy*(f) 

Gyx(f> 

Gxx(f> 

Gyy(f) 

H(f ) 

«yx(f> 

!Kt) 

K 

N 

HK 


amplitude  at  mesh  frequency  and  harmonics 
original  signal  after  time  synchronous  averaging 
analytic  signal 
difference  signal 

discrete  value  In  the  difference  signal,  D(t) 
mean  value  of  the  difference  signal,  D<t) 
discrete  value  in  a  frequency  array 
discrete  value  In  a  time  record 

discrete  frequency  array  (discrete  fourler  transform  of  a  tire 
record) 

discrete  time  record 

frequency  domain  function  (fourler  transform  of  a  time  domain  func¬ 
tion) 

frequency  (cycles/second) 

minimum  resolution  frequency 

first  order  sidebands  level 

cross  power  spectrum 

complex  conjugate  of  Gxy<f> 

cross  power  spectrum 

auto  power  spectrum  of  Input 

auto  power  spectrum  of  output 

frequency  response  function 

estimate  of  frequency  response  function 

unit  Impulse  response  function 

kurtosls,  fourth  moment  of  an  array  of  values 

number  of  data  points  in  a  discrete  time  record 

normalized  kurtosls 
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NS 


number  of  samples 
Nt  number  of  time  averages 

PL  peak  level  of  signal 

PP  peak  to  peak  level  of  signal  average 

R(t)  regular  meshing  components  of  signal  average 
RMS  standard  deviation  of  signal  average 

RMSDS  standard  deviation  of  difference  signal,  D(t) 

RMSRC  standard  deviation  of  regular  meshing  components  of  signal,  R(t> 

S/N  signal  to  noise  ratio 

T  length  of  time  record  (seconds) 

t  time  (seconds) 

w  frequency  (radians/second) 

X(f )  fourier  transform  of  Input  signal 

X*(f)  complex  conjugate  of  X(f) 

Y(f)  fourler  transform  of  output  signal 
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